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Summary. We present a comprehensive review of the physics of hadron and jet production 
at large transverse momentum in high-energy nucleus-nucleus collisions. Emphasis is put on 
experimental and theoretical "jet quenching" observables that provide direct information on 
the (thermo)dynamical properties of hot and dense QCD matter. 



1 Introduction 

The research programme of high-energy nucleus-nucleus physics is focused on the 
study of the fundamental theory of the strong interaction - Quantum Chromo Dy- 
namics (QCD) - in extreme conditions of temperature, density and small parton 
momentum fraction (low-x) - see e.g. [11 for a recent review. By colliding two 
heavy nuclei at ultrarelativistic energies one expects to form a hot and dense de- 
confined medium whose collective (colour) dynamics can be studied experimentally. 
Lattice QCD calculations |2 1 predict a new form of matter at energy densities (well) 
above e„„ « 1 GeV/fm^^ consisting of an extended volume of deconfined and chirally- 
symmetric (bare-mass) quarks and gluons: the Quark Gluon Plasma (QGP) |3 1. 

Direct information on the thermodynamical properties (like temperature, energy 
or particle densities, ...) and transport properties (such as viscosity, diffusivity and 
conductivity coefficients) of the QGP can be obtained by comparing the results for 
a given observable <Paa measured in nucleus-nucleus (AA, "QCD medium") to those 
measured in proton-proton (pp, "QCD vacuum") collisions as a function of centre- 
of-mass (cm.) energy y/s^, transverse momentum pj, rapidity y, reaction centrality 
(impact parameter b), and particle type (mass m). Schematically: 

„ , , — "hot/dense QCD medium" (t>AAi \/s^,PT ,y,m;b) 

RAA{V^,PT,y,m;b) ^ — — °^ —f—Tr ^ 

QCD vacuum 9pp{^ s ,pT ,y,m) 

Any observed enhancement {Raa > 1) and/or suppression {Raa < 1) in this ratio can 
then be linked to the properties of strongly interacting matter after accounting for 
a realistic (hydrodynamical) modeling of the space-time evolution of the expanding 
system of quarks and gluons produced in the collision. 

This article presents an up-to-date review of experimental and theoretical studies 
of quantities such as Raa{ ^/s^,PT,y,iii;b) forhigh-pr hadrons and jets produced at 
RHIC collider energies and prospects for the LHC. 
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2 Jet quenching and parton energy loss in QCD matter 

In this first Section, we introduce the general concepts, basic variables and formulas 
of energy loss of a fast charged particle in a dense thermalised plasma (starting with 
the somewhat simpler QED case), and we enumerate the expected phenomenological 
consequences of QCD energy loss for gluons and light- and heavy-quarks traversing 
a hot and dense QGP. 

2.1 Hard probes of hot and dense QCD matter 

Among all available observables in high-energy nuclear collisions, particles with 
large transverse momentum and/or mass, pT,in ^ Qq S> A^^.^, where Qq = ^(1 GeV) 
and Ag^^ « 0.2 GeV is the QCD scale, constitute valuable tools to study "tomo- 
graphically" the hottest and densest phases of the reaction (Fig. [TJ. Indeed, such 
"hard probes"(i) originate from partonic scatterings with large momentum transfer 
and thus are directly coupled to the fundamental QCD degrees of freedom, (ii) are 
produced in very short time-scales, x ^ l/pr l/Qg ^ 0.1 fm/c, allowing them to 
propagate through (and be potentially affected by) the medium, and (iii) their cross 
sections can be theoretically predicted using the perturbative QCD (pQCD) frame- 
work im. 



Fig. 1. Examples of hard probes whose modifications in high-energy AA collisions provide 
direct information on properties of QCD matter such as the transport coefficient q, the initial 
gluon rapidity density dN^ /dy, and the critical temperature r„i, and energy density e„„ yj. 

Jet production in hadronic collisions is a paradigmatic hard QCD process. An 
elastic (2 — > 2) or inelastic (2 — > 2 + X) scattering of two partons from each one of 
the colliding hadrons (or nuclei) results in the production of two or more partons 
in the final-state. At high pr, the outgoing partons have a large virtuality Q which 
they reduce by subsequently radiating gluons and/or splitting into quark-antiquark 
pairs. Such a parton branching evolution is governed by the QCD "radiation prob- 
abilities" given by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions down to virtuaUties l GeV^). At this point, the produced partons frag- 
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ment non-perturbatively into a set of final-state hadrons. The characteristic colli- 
mated spray of hadrons resulting from the fragmentation of an outgoing parton is 
called a "jet". 




(quenched] jet 



Fig. 2. "Jet quenching" in a head-on nucleus-nucleus collision. Two quarks suffer a hard scat- 
tering: one goes out directly to the vacuum, radiates a few gluons and hadronises, the other 
goes through the dense plasma created (characterised by transport coefficient q, gluon density 
dN^ /dy and temperature T), suffers energy loss due to medium-induced gluonstrahlung and 
finally fragments outside into a (quenched) jet. 



One of the first proposed "smoking guns" of QGP formation was "jet quench- 
ing" f6l i.e. the attenuation or disappearance of the spray of hadrons resulting from 
the fragmentation of a parton having suffered energy loss in the dense plasma pro- 
duced in the reaction (Fig. |2|i. The energy lost by a particle in a medium, AS, pro- 
vides fundamental information on its properties. In a general way, AE depends both 
on the characteristics of the particle traversing it (energy E, mass ;«, and charge) and 
on the plasma properties (temperature T, particle-medium interaction coupling a, 
and thickness L), i.e. AE{E,m, T,a,L). The following (closely related) variables are 
extremely useful to characterise the interactions of a particle inside a medium: 

• the mean free path A- = 1 / (per), where p is the medium density (p °^ for an 
ideal gas) and o the integrated cross section of the particle-medium interactior0, 

• the opacity N — L/X or number of scatterings experienced by the particle in a 
medium of thickness L, 

• the Debye mass moiT) ^ gT (where g is the coupling parameter) is the inverse of 
the screening length of the (chromo)electric fields in the plasma, mo characterises 
the typical momentum exchanges with the medium and also gives the order of 
the "thermal masses" of the plasma constituents, 

• the transport coefficient q = m^/X encodes the "scattering power" of the medium 
through the average transverse momentum squared transferred to the traversing 
particle per unit path-length, q combines both thermodynamical (mo, p) and dy- 
namical (a) properties of the medium ||7]|8]|9l: 

q = ml[)/X = mj) pa. (2) 



' The QED and QCD coupling "constants" are aem = I (43t) and a, = g^/(47l) respectively. 
^ One has A. ~ (ar)^' since the QED, QCD screened Coulomb scatterings are Q^i (x/T^ . 
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As a numerical QCD exampl^E let us consider an equilibrated gluon plasma 
at r = 0.4 GeV and a strong coupling « 0.5 jlOl . At this temperature, the 
particle (energy) density is p^. ^ 16/%^ C(3) -T^ ~ 15 fm"^ (e^ = S%^/15 ■ 
w 17 GeV/fm^), i.e. 100 times denser than normal nuclear matter (p = 0.15 fm^^). 
At leading order (LO), the Debye mass is mo = (47tai)'/^r p^i 1 GeV. The LO 
gluon-gluon cross section is <5gg ~ 97ta^/(2OT^) « 1.5 mb. The gluon mean free 
path in such a medium is A,g = 1 /{^g<3gg) ~ 0.45 fm (the quark mean-free-path is 
'^q = CA/CF'kg~ 1 fm, where Ca/Cf = 9/4 is the ratio of gluon-to-quark colour 
factors). The transport coefficient is therefore q ~ m^/Xg ~ 2.2 GeV^/fm. Note 
that such a numerical value has been obtained with a LO expression in a, for 
the parton-medium cross section. Higher-order scatterings (often encoded in a 
"/T-factor" « 2 - 4) could well result in much larger values of q. 
• the diffusion constant D, characterising the dynamics of heavy non-relativistic 
particles (mass M and speed v) traversing the plasma, is connected, via the Ein- 
stein relations 

D - 2rVK =r/(M rio) (3) 

to the momentum diffusion coefficient K - the average momentum squared gained 
by the particle per unit-time (related to the transport coefficient as K « ^ v) - and 
the momentum drag coefficient r[o- 

2.2 Mechanisms of in-medium energy loss 

In a general way, the total energy loss of a particle traversing a medium is the sum of 
collisional and radiative termO ~ AEcoii +^E,-aci- Depending on the kinematic 
region, a (colour) charge can lose energ}0 in a plasma with temperature T mainly by 
two mechanismo 




Fig. 3. Diagrams for collisional (left) and radiative (right) energy losses of a quark of energy 
E traversing a quark-gluon medium. 



• Collisional energy loss through elastic scatterings with the medium constituents 
(Fig. |3] left) dominates at low particle momentum. The average energy loss in 

^ For unit conversion, multiply by powers of fic ~ 0.2 GeV fm (other useful equalities: 
10 mb = 1 fm^ and 1 GeV^^ = 0.389 mb). 

^ In addition, synchrotron-, Cerenkov- and transition-radiation energy losses can take place 
respectively if the particle interacts with the medium magnetic field, if its velocity is greater 
than the local phase velocity of light, or if it crosses suddenly from one medium to another. 
Also, plasma instabilities may lead to energy losses. Yet, those effects - studied e.g. in 111! 
I12II13I I141 for QCD plasmas - are generally less important in terms of the amount of Eio^s- 

^ Note that if the energy of the particle is similar to the plasma temperature, E ~ ff{T), the 
particle can also gain energy while traversing it. 

* Note that the separation is not so clear-cut since the diagrams assume well-defined asymp- 
totic out states, but the outgoing particles may still be in the medium and further rescatter. 
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one scattering (with cross section do/dt, where t — is the momentum transfer 
squared) in a medium of temperature T, is: 

(AE^^if)^— t—dt. (4) 

\ ' or X„2 dt 

• Radiative energy loss through inelastic scatterings within the medium (Fig. [31 
right), dominates at higher momenta. This loss can be determined from the cor- 
responding single- or double-differential photon or gluon Bremsstrahlung spec- 
trum ((0 diraci / d(x) or 03d^Irad/d(odk^, where (O, are respectively the energy 
and transverse momentum of the radiated photon or gluon): 

<"=/'«;^*, 0, A.;.r=/7"-«ie^*.4.(5) 

For incoherent scatterings one has simply: AE'"' — • AE^"'^"', where = L/X is the 
medium opacity. The energy loss per unit length or stopping powe^is: 



dE _ {AE" 



dl L 

which for incoherent scatterings reduces to: —dE/dl = ^A/i'*™') /A,. 



(6) 



Energy losses in QED 

As an illustrative example, we show in Fig. |4]the stopping power of muons in cop- 
per. At low and high energies, the collisional (aka "Bethe-Bloch") and the radiative 
energy losses dominate respectively. 
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Fig. 4. Stopping power, —dE/dl, for positive muons in copper as a function of |3y= p/Mc (or 
momentum p). The solid curve indicates the total stopping power [15 1. 



Yet, the hot and dense plasma environment that one encounters in "jet quench- 
ing" scenarios is not directly comparable to the QED energy loss in cold matter 
represented in Fig.|4] A recent review by Peigne and Smilga lfT6l presents the para- 
metric dependences of the energy loss of a lepton traversing a hot QED plasma with 

^ By 'stopping power', one means a property of the matter, while 'energy loss per unit length' 
describes what happens to the particle. For a given particle, the numerical value and units 
are identical (and both are usually written with a minus sign in front). 
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temperature T and Debye-mass mo. In a simplified manner, inserting the Coulomb 
(lepton-lepton) and Compton (lepton-photon) scattering cross sections in Eq. (|4|i and 
using Eq. (|6]l, one obtains the collisional energy losses per unit-length: 

• Light lepton (M^ < ET): ^o?T^ In (^^) - a m^, In (^^) 
. Heavy lepton (M^ » ET): -'-Mml ^ 2^^2j^2 i„ _ „ ,„2 i„ 

For radiative losses, the amount of photon emission, depends chiefly on the thickness 
of the plasmc0. For thin media (L <C the traversing particle suffers at most one 
single scattering and the QED radiation spectrum is just given by the Bethe-Heitler 
(BH) bremsstrahlung expression ifTTI . On the contrary, for thick media (L ^ A,) there 
are A^(=opacity) scatterings and the Landau-Pomeranchuk-Migdal (LPM) 1 1 8 1 coher- 
ence effect reduces the amount of radiation compared to times the BH spectrum. 
Making use of Eq. Q one obtains, from each corresponding radiation spectra, the 
following parametric expressions lfT6l : 

• "Bethe-Heitler'0(i< A,): ^ a (L^ml/X) ■ (o/E^ a qL^ ■ &/E^ 

AEZ-^qL^-o?T'L^ =^ = ^ ^ a' t' L (7) 

• Landau-Pomeranchuk-Migdal (L > A.) : &'-^^a^L Jt^ In {E^/{&T)) 



AeZ"" --a^L^ETHn{E/T) ^ -^^^a^ ^ ET^ HE/T) (8) 

We see that, in general, the radiative energy losses of an energetic lepton crossing 
a hot QED plasma are much larger than their collisional losses. Yet, if the particle 
is heavy, the amount of radiation at angles within a cone Q < M/E is suppressed 
by a factor mj)/M^ ("dead cone" effect, see later) resulting in a reduction of the 
bremsstrahlung emission by a factor m^/hfi ^ aT^/M^. 



Energy losses in QCD 

The main differences between QED and QCD energy losses result from the non- 
Abelian nature of QCD: i.e. the fact that gluons can also interact with themselves (at 
variance with photons in QED) introduces several important changes. First, the QCD 
coupling a.s runs more rapidly than (Xem (at least for not asymptotically-high temper- 
atures), and the scale Q at which {Q) is evaluated needs to be explicitly considered 
in all calculations of collisional energy losses fT9^, "201. Second, it is crucial to take 
into account the different coupling of quarks and gluons with the medium. The rela- 
tive strengths of the three distinct QCD vertices: (XsCp for q ^ qg, UsCa for g ^ gg, 
and UsTf for g ^ qq, are determined by the structure (Casimir factors Cr) of the 
gauge group describing the strong forcg^H The probability for a gluon (quark) to ra- 
diate a gluon is proportional to the colour factor d = 3 (C/r = 4/3). In the asymptotic 

^ We consider here the formulas where the charged particle is produced inside the plasma, as 
this is the typical situation encountered in a QGP. 

' The LPM effect describes the fact that, since it takes a finite time to emit a photon, neigh- 
bouring medium particles interfere coherently (destructively) and act as one effective scat- 
tering centre, inducing single photon radiation. 

Strictly speaking, as discussed in 1 16], for light particles produced inside a plasma, there 
is no BH regime whatsoever: the L'-dependence appearing in AE*^ is due to in-medium 
formation-time constraints absent in the truly (asymptotic) BH regime. 
" Foi SU{Nc) wilhNc the number of colours: Ca =Nc, Cp = {N^-l)/2Nc and Tf = 1/2 (21]. 



Jet quenching 7 



limit, in which the radiated gluons carry a small fraction of the original parton mo- 
mentum (and neglecting the gluon splitting into quark-antiquark pairs, proportional 
to the smaller factor = 111), the average number of gluons radiated by a gluon is 
Ca/Cf = 91 A times higher than that radiated by a quark (that is the reason why gluon 
jets have a larger (and softer) hadron multiplicity than quark jets). 

QCD collisional energy loss 

The collisional energy loss due to elastic scattering of a parton of energy E in- 
side a QGP of temperature T was originally estimated by Bjorken f6l and Braaten- 
Thoma ll22l and later improved (including running coupling, finite energy kinemat- 
ics, and quark-mass effects) by various authors ||T9l |20l |23l . Using Eq. (|4]i with 
the momentum-transfer integral limits given by (i) the QGP Debye-mass squared 
tmin = ^{T) ^ A%asT^{l +Nf/6) and (ii) t^ax ~s ^ ET, and taking the dominant 
contribution to the parton-parton /-differential elastic cross section 

— « Ci , with av(f) = ^ — —-5 — - (9) 

dt ' f2 ' (33-2«/)ln(f/A2gc^) 

where C, = 9/4, 1,4/9 are the colour factors for gg, gq and qq scatterings respec- 
tively, one finally obtains (forE > /T) |20|: 

• Light-quark, gluon: ^ = \CRa,{ET) ml In {^■ 

• Heavy-quai-k:-^%^L = -^%^| - 2 ^^^^ ^,^2^„ /j,t^m„ / irr 



a,{M^)a,{ET) In 



dl dl \ci 9' 

with Cr = 4/3 (3) being the quark (gluon) colour charge. The amount of AEcuii is 
linear with the medium thickness, and it depends only logarithmically on the initial 
parton energy. As a numerical example, taking E = 2Q GeV, M = 1 .3 GeV/c^ (charm 
quark) and a medium with T — QA GeV and nio = 1 GeV/c^, the elastic energy losses 
per unit-length are —dEcii /dl\^ — 2.3 GeV/fm and —dEcon/dl\^ = 2.6 GeV/fm. 



QCD radiative energy loss 



i'q^qg(z) - Cf ^ ^ , and Pg^gg(z) = Q ' _Z ' dO) 



The dominant mechanism of energy loss of a fast parton in a QCD environment is 
of radiative nature ("gluonstrahlung") 1241 l25l l26l l27l : a parton traversing a QGP 
loses energy mainly by medium-induced multiple gluon emission. The starting point 
to determine the radiation probabilities in QCD are the DGLAP splitting functions 
in the vacuum |5 1 

z(l-z) 

(where z = (o/E is the fraction of energy of the parent parton taken by the radiated 
gluon), modified to take into account the enhanced medium-induced radiation. The 
resulting radiated gluon spectrum, co dlrad /doi Pq i,^^(^(o/ E^ , has been computed 
by various groups under various approximations (see Section[32|. All medium mod- 
ifications are often encoded into the "transport coefficient" parameter, q, introduced 
previously, see Eq. (|2]i. For thin (thick) media, i.e. for L-^X(L^X), one deals with 
the Bethe-Heitler (Landau-Pomeranchuk-Migdal) gluonsstrahlung spectrum. In the 
LPM case, one further differentiates between the soft or hard gluon emission cases 
with respect to the characteristic gluonstrahlung energ}{3 ra^ = ^q L^. Making use 
of Eq. (|5]), the basic QCD radiative energy loss formulas read [ 16J: 

• "Bethe-Heitler'[3 (BH) regime (L < A.): 

'■^ Up to prefactors, cOc is the average energy lost in the medium (for cOe < E): cOc ~ 
2(AE,.ad) /{asCn). 

'■^ In reahty, BH usually refers to the single scattering spectrum of an asymptotic particle 
produced outside the medium (see footnote for the analogous Bethe-Heitler QED formula). 
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co^»a,^LVco =^ AEf^^^,^asqLHn{E/{mlL)) (11) 
Landau-Pomeranchuk-Migdal (LPM) regime (L 3> A,): 



We note two things. First, because of the destructive interference, the LPM spectrum, 
^dlrad /d(0 °= CO^ 'Z^, is Suppressed in the infrared (i.e. for small co's) compared to the 
independent Bethe-Heitler gluon spectrum, (Hdlrad / d(x) °^ (x)^K Note also that, due to 
the steeply falling spectrum of the radiated gluons, the integrated LPM energy loss is 
dominated by the region co < co,. Second, the QCD energy loss shows a characteristic 

dependence on the plasma thickness which however, as noted in [161, is also 
present in the case of Abelian (QED) plasmas, see Eq. O, and is a general feature 
of the medium-induced energy loss of any in-medium newborn particle. The main 
distinctions of the energy loss in a QCD compared to a QED plasma are the presence 
of different colour factors for the q and g charges (AErad °^q°^ '^panicle-medium which 
in the QCD case is proportional to Cr) and the extra logarithmic dependence of ls.E,-ad 
on the energy E of the traversing particle. 

For a gluon with E = 20 GeV in a medium with q = 2 GeV^/fm and L = 6 fm, 
dErad/dl is ^(10 GeV/fm) (to be compared with the elastic losses of ff{2 GeV/fm) 
estimated before). As seen in Fig. |5] for a more realistic phenomenological case, 
Afico/; is in general a small correction compared to l^rad for light quarks and gluons 
but it can be an important contribution for slower heavy-quarks (see next). 




Fig. 5. Comparison of the average radiative and elastic energy losses of light-quarks (left) 
and light- and heavy-quarks (right) passing through the medium produced in central AuAu 
colUsions at RHIC energies as obtained by the AMY (28| and DGLV (29) models (see later). 



Heavy-quark radiative energy loss ( "dead cone " effect) 

Gluon bremsstrahlung off a heavy quark differs from that of a massless parton, al- 
ready in the vacuum. Due to kinematics constraints, the radiation is suppressed at 
angles smaller than the ratio of the quark mass M to its energy E. The double- 
differential distribution of gluons of transverse momentum k± and energy CO ra- 
diated by a heavy quark at small-angles {k± « CO 9), differs from the standard 
bremsstrahlung spectrum by the factor 

dIrad.Q ^sCf k\ dl,.ad f , , ^lY^ „ M 1 
CO ^- = —5 \, T,. W CO T- • l + 'TT ,60 = — = -. (13) 

d(odkl n (yt2 +0)292)2 dcodk^ \ J E y 
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This effect, known as the "dead cone" |30|, results in a reduction of the total gluon 
radiation emitted off heavy-quarks. In the medium, the total amount of reduction 
depends on a non-trivial way on the various scales {E,M,L) of the problem [ 16 1. In a 
simplified way, the effect is 6'{mo/M) - compared to ff{m^/M^) in the QED case - 
i.e. for a plasma with Debye-mass nio = 1 GeV/c^, the reduction of radiative energy 
loss for a charm (bottom) quark of mass 1.3 (4.2) GeV/c^ is of order -25% (75%). 



2.3 Phenomenological consequences of parton energy loss 

Medium-induced parton energy loss in AA reactions results in various observable 
consequences compared to the same "free space" measurements in proton-proton 
ipp) collisions. The presence of jet quenching manifests itself via : 



(i) a suppression of the spectrum (dNAA/dpr) of high-y^j- hadrons | 

(ii) unbalancedback-to-backhigh-/??- dihadron azimuthal correlations {dNpair/d^) OTl . 

(iii) modified energy-particle flow (softer hadron spectra, larger multiplicity, in- 
creased angular broadening, ...) within the final jets ||32l [33l l34l |35l [36l [37l . 

In addition, due to the aforementioned hierarchy of flavour-dependent radiative en- 
ergy losses: AEnuiig) > ^Erad{q) > ^Erad{c) > AErad{b), all these medium effects 
are expected to be larger for gluons and u, d, s quarks than for c orb quark0. 

(i) High-pT (leading) hadron spectra: 

The leading hadron of a jet is the hadron that carries the largest fraction of the energy 
of the fragmenting partorF^ In a heavy-ion collision, if the parent parton suffers 
energy loss, the energy available for such hadrons is reduced and consequently their 
spectrum is depleted compared to pp. From the measured suppression factor one can 
determine A£/o„ (see Eq. ( [36] l later) and estimate various properties of the produced 
plasma such as: 

the average transport coefficient {q), from Eqs. (fTTT i. ( fT2] i via {AE) ^ {q)L^, 



• the initial gluon density dN^ /dy of the expanding plasma (with original trans- 



verse area Ax = T^R^ ~ 150 fm^ and thickness L), from 1261 : 



AEc<a^CR^^L. (14) 

Aj_ dy 



High-pT spectra in pp and AA collisions are discussed in detail in Chapter |4ni 
(ii) High-pT di-hadron correlations: 

Parton-parton 2^2 scatterings are balanced in pj i.e. they are back-to-back in 
azimuthal angle (A(|) « n). Such azimuthal correlation is smeared out if one or 
both partons suffer rescatterings ("kr broadening") in a dense plasma. The dijet- 
acoplanarity arising from the interactions of a parton in an expanding QGP is 
i^T)med ~ i^%/'^)E^^{L/iQi) °^ qL |[39l and, thus, the final azimuthal correlations 
between the hadrons issuing from quenched partons will show a dependence on the 
transport coefficient and thickness of the medium: d^Npair/dA(^ = f{q,L). 

In addition, it has been proposed that a parton propagating through a QGP with 
supersonic (p > c,) or "superluminal" (P > 1 /n) velocities can generate a wake of 
lower energy gluons with either Mach- ||40] 111] |42l or Cerenkov-like ||42l |43] gU 
conical angular patterns. Such a conical emission can propagate, after hadronisation, 
into the final correlations of the measured hadrons with respect to the jet axis: 



The heaviest top-quark decays into W b immediately (x < 0. 1 fm/c) after production. 
High-pj- inclusive hadron spectra are dominated by particles with (z) ~ 0.5 - 0.7 1381 
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• In the first case, the speed of sound of the traversed matteiF^ — dP /d^, can 
be determined from the characteristic Mach angle 9^ of the secondary hadrons: 

cos(0m) = , (15) 

• In the second scenario, the refractive index of the medium, n « ^Jz'r where Er 
is the gluon dielectric constant, can be estimated from the Cerenkov angle of 
emission 9^ of the hadrons: 

cos(9,) = \ = (16) 

High-/??- hadron correlations in nucleus-nucleus collisions are covered in Chapter|5] 

(Hi) Jet spectra, jet shapes, and fragmentation functions: 

The measurements of fully reconstructed (di)jets or of jets tagged by an away-side 
photon ||451 or Z-boson |46 | in heavy-ion collisions allow one to investigate - in 
much more detail than using single- or double- hadron observables - the mechanisms 
of in-medium parton radiation as well as to obtain the transport coefficient q via: 

• Medium-modified jet profiles ll32l l47l . through the differential f)'""'{r;q) and 
integrated '^"^'^{r;q) jet-shapes (r is the distance to the jet axis) and the thrust 
.^{q) variables, 

• Medium-modified fragmentation functions 1481. ^I^arton^hadmn 

{z,Q^), where 

z = P hadron I P parton is the fraction of the jet energy carried by a hadron, which 
very schematically can be written as 

^TM{z-A,L) c P{z-q.L) ® 9J^l{z) , (17) 

where the correction P{z,q,L) can be connected to the QCD splitting functions, 
Eq. ( [Tol l, modified according to various possible prescriptions (see Section [33T l 
to take into account medium-induced gluon radiation. 

Jet physics studies in heavy-ion collisions are covered in Chapter|6] 



3 Parton energy loss phenomenology 

The use of fast partons as calibrated tomographic probes of hot and dense QCD mat- 
ter in heavy-ion collisions relies on the possibility to compute theoretically (i) their 
perturbative production cross sections, and (ii) their modifications suffered while 
propagating through a strongly-interacting medium. We discuss here the basic pQCD 
principles used to compute high-pr hadron (and jet) cross sections, and then we out- 
line the various existing parton energy loss schemes. 

3.1 QCD factorisation in high-/9r hadron and jet production in AA colUsions 

Because of asymptotic freedom, the QCD coupling is small for high-energy (short 
distance) parton interactions: a.s (G^ —^oo^^O. The single inclusivg ^1 production of 
a high-pT parton c in a parton-parton collision, ab ^ c +X, can thus be computed 

The speed-of-sound - namely the speed of a small disturbance through the medium - for 
an ideal QGP (with e = 3P, where P is the pressure) is simply = 1 / \/3. 
Inclusive refers to the consideration of all possible channels that result in the production of 
a given particle c, without any particular selection of the final-state X. 
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using perturbation theory techniques. Over short distances, the infinite number of 
Feynman diagrams that would theoretically result in the production of the outgoing 
parton c, can be approximated accurately by a much more manageable number of 
terms. In high-energy hadwn-hadwn collisions, the production of high-pj particles 
can be computed from the underlying parton-parton processes using the QCD "fac- 
torisation theorem" (see sketch in Fig. |6]l [491. The production cross section of a 
h\g\\-pT hadron h can be written as the product 

where (Sab^cxi^i^^i^Q^) is the perturbative partonic cross section computable up to 
a given order in tts, and where the two non-perturbative terms: 

• fa/ A {^tQ^)'- parton distribution functions (PDF), encoding the probability of find- 
ing a parton of flavour a and momentum fraction x = p parton/ Pmicieus inside the 
nucleus A, 

• S>c^h{zTQ^)'- fragmentation function (FF), describing the "probability" that the 
outgoing parton c fragments into a final hadron h with fractional momentum 

Z — P hadron I P part 0117 

are universal (i.e. process-independent) objects that can be determined experimen- 
tally e.g. in deep-inelastic e*-nucleus and e+e^ collisions, respectively. For the case 
of total parton (i.e. jet) cross section, one simply sets S>c^h = 8(1 — z) in Eq. (fTsT l. 




hadrons 

Fig. 6. Sketch of dijet production and pQCD coUinear factorisation in hadronic collisions: 
fa/ A (-'■) the PDFs, S>i^k (z) the FFs, and ISR (FSR) represents initial (final)-state radiation. 

The basic assumption underlying the factorised form of Eq. ( fTsT l is that the char- 
acteristic time of the parton-parton interaction is much shorter than any long-distance 
interaction occurring before (among partons belonging to the same PDF) or after 
(during the evolution of the struck partons into their hadronic final-state) the hard 
collision itself. The validity of Eq. ( fTsT l holds thus on the possibility to separate long- 
and short-distance effects with independent QCD time- (length-) scales, as well as on 
the "leading-twist'[3 assumption of /nco/ierenf parton-parton scatterings. Since par- 
tons are effectively "frozen" during the hard scattering, one can treat each nucleus as 
a collection of free partons. Thus, with regard to high pj production, the density of 



Processes in which more than one parton from the same hadron/nucleus interact coherently, 
are called "higher-twist" processes. 
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partons in a nucleus with mass number A is expected to be simply equivalent to that 
of a superposition of A independent nucleona'l: fa/Ai^^Q^) ~A ■/«/a'(^i2^)- Thus, 

do'/^ih~A-B ■faip{x,Q^)® fblp{x.Q^)<i^dc'jic<^%^h{z.Q^). (19) 

From (fTSl l, it is clear that QCD factorisation implies that hard inclusive cross sections 
in a A B reaction scale simply as A • B times the corresponding pp cross sections: 

dcs%'^ =A-B-do^^f. (20) 

Since nucleus-nucleus experiments usually measure invariant yields for a given cen- 
trality bin (or impact parameter b), one writes instead: 

dN^T'ib)^{TAB{b))-da';^';'', (21) 

where the nuclear overlap function at b, TAsib), is determined within a geometric 
Glauber eikonaj^ model from the measured Woods-Saxon distribution for the in- 
teracting nuclei (51). Intuitively, one can think of the nuclear-overlap 7aa(^) as a 
function that characterises the surface profile of two "beams" of nucleons colliding 
at a distance b. The units of [area]^' of Taa indicate that it represents somehow 
the effective "parton (integrated) luminosity" of the collision. Since the number of 
inelastic nucleon-nucleon (A^A^) collisions at b, Ncoii{b), is proportional to TAsib): 
Ncoii{b) = TAB{b) ■ o"Jf', one also writes often Eq. (l2Tl i as: 

dN'A'S'ib) = (A^,„„ {b)) ■ dN'^f. (22) 

For minimum-bia A B collisions, the average nuclear-overlap and number of A^A^ 
collisions take a simple forrr0: (Tab) —^B/af^ and (Ncuii) =AB- Onn/o^ab- The 
standard method to quantify the effects of the medium on the yield of a hard probe 
in a AA reaction is thus given by the nuclear modification factor. 

D i u\ d^NAA/dydpT 

RAA{PT,y;b) = , (23) 

{TAA(b)) X d'^Opp/dydpT 

This factor - a quantitative version of the ratio sketched in Eq. ([T]i - measures the 
deviation of AA at b from an incoherent superposition of A^A^ collisions (Raa = !)■ 
This normalisation is often known as "binary collision scaling". 



3.2 Jet quenching models 

The energy-loss formulas presented in Section 12.21 refer to an idealistic situation 
with an infinite-energy parton traversing a static and uniform QGP with an ideal- 
gas equation-of-state (EoS). Experimentally, the situation that one encounters with 
realistic plasmas in heavy-ion collisions is more complex: 

• first, there is no direct measurement of the traversing parton but (in the best case) 
only of the final-state hadrons issuing from its fragmentation, 

• the traversing partons can be produced at any initial point within the fireball and 
their energy spectrum is steeply (power-law) falling. 

In reality, nuclear PDFs are modified compared to proton PDFs by initial-state 

"(anti)shadowing" effects (see 1501 for a recent review). 
^'^ The 'eikonal' approximation assumes that the particle trajectories are simple straight lines. 
^' Minimum-bias collisions are those where there is no specific selection of the final-state (e.g. 

in particular for heavy-ions, no centrality selection). 
22 E.g. fotAuAu at /i;^ = 200 GeV (<''' = 41 mb, of °„ = 7000 mb): {TauAu) = 5.5 (23.3) 

mb^' and (Ncoii) = 230 (955) for minimum-bias (10% most central) collisions, resp. 
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• the temperature and density of the plasma, and correspondingly its Debye-mass 
and transport coefficient, are position-dependent: mo(r),^(r), 

• the produced plasma is expanding with large longitudinal (transversal) velocities, 
P w 1 (0.7), and so the medium properties are also time-dependent: m[){l),q{l), 

• the finite-size of the medium and associated energy loss fiuctuations, have to be 
taken into account. 

All those effects can result in potentially significant deviations from the analytical 
formulas of Section IZ21 (e.g. in an expanding plasma the dependence of AEi-ad on 
the medium thickness L, Eq. ( fTSl i. becomes effectively linear, Eq. ( fT4b . rather than 
quadratic). Four major phenomenological approaches (identified often with the ini- 
tials of their authors) have been developed ll52l to connect the QCD energy loss 
calculations with the experimental observables mentioned in Section l23] 

• Path-integral approach to the opacity expansion (BDMPS-LCPI/ASW) |l53] l54l 
|55]|56l|25l|521|27l|58l|59l 

• Reaction Operator approach to the opacity expansion (DGLV) |l60l |26] [61] |62] |29l 

• Higher Twist (HT) 1 63 64 65] |66l El ISH) 

• Finite temperature field theory approach (AMY) 1691170117111721 

The models differ in their assumptions about the relationships between the relevant 
scales (parton energy E and virtuality Q^, and typical momentum /u ~ and spatial 
extent L of the medium), as well as by how they treat or approximate the space-time 
profile of the medium. In practical terms, all schemes are based on a pQCD factorised 
approach, i.e. on Eq. (ITSl l. where the entire effect of energy loss is concentrated 
on the calculation of the medium-modified parton fragmentation functions into final 
hadrons: (z) ^c™*/' {z' , q) ■ The final hadronisation of the hard parton is always 
assumed to occur in the vacuum after the parton, with degraded energy (z' < z), has 
escaped from the system (Fig.|7]l. 



Fig. 7. Schematic representation of parton energy loss in a QGP, implemented via rescaling of 
the energy of the traversing parton at the point where it fragments into hadrons 1481 . 



BDMPS-LCPI & ASW 

The approaches of Baier, Dokshitzer, Mueller, Peigne and Schiff (BDMPS) 1561 |25l 
|73]|74l and the Light-Cone-Path-Integral (LCPI) by Zakharov |53 | compute energy 
loss in a coloured medium in a multiple soft-scatterings approximation. A hard par- 
ton traversing the medium interacts with various scattering centres and splits into an 
outgoing parton as well as a radiated gluon (Fig.|8]l. The propagation of the traversing 
parton and radiated gluons is expressed using Green's functions which are obtained 
by a path integral over the fields. The final outcome of the approach is a complex an- 
alytical expression for the radiated gluon energy distribution (odI/d(x) as a function 
of the transport coefficient q, Eq. (|2]i, defined perturbatively as llTSl : 



where p is the medium density of scattering centres (mainly gluons), k± is the trans- 
verse momentum of the radiated gluon, and da is the differential parton-medium 




h 




(24) 
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E parton 




Fig. 8. Typical gluon radiation diagram in tiie BDMPS approach! |25|. 



cross section. The medium-modified parton-to-hadron fragmentation functions are 
modelled as 

(25) 



^med 



where the quenching weights q) - computed by Armesto, Salgado and Wiede- 
mann (ASW) ||59l l76l IJTI - encode the probability (assumed Poissonian) that the 
propagating parton loses a fraction of energy e = AE/E due to n gluon emissions 



n „ 

n / d(Oi 



dr""'{q) 



5 e 




exp 



d(£)- 



dl 



med 



dm 



(26) 



The quenching weights have been implemented in a Monte Carlo model, the Parton- 
Quenching-Model (PQM) |78t |79„l accounting for a realistic description of the parton 
production points in a static QGP. The transport coefficient q is used as the fit param- 
eter for the data. The longitudinal expansion of the plasma can be taken into account 
by rescaling the transport coefficient according to the following law [58 1: 



2 /"^o 



di (t-to) q{i) 



(27) 



where q{'i) — q{xo) {xo/z)"" and a characterises the time-dependence of the plasma 
density: p(t) o= t^". A purely longitudinal (or Bjorken) expansion corresponds to 
a = 1, and is often assumed in phenomenological applications. When To ^ L, 
Eq. dzTl i reduces to (q) ~ ^(to)2to/L | 



DGLV 



The Gyulassy-Levai-Vitev (GLV) iOl |60l EH IMl (aka DGLV Mi^) approach 
calculates the parton energy loss in a dense deconfined medium consisting, as in the 
BDMPS approach, of almost static (i.e. heavy) scattering centres (Fig. |9]l produc- 
ing a screened Coulomb (Yukawa) potential. Although both approaches are equiv- 
alent ||27]|59l|82l, at variance with the BDMPS mM/fi>Ze-so/f bremsstrahlung, GLV 
starts from the single-hard radiation spectrum which is then expanded to account 
for gluon emission from multiple scatterings via a recursive diagrammatic proce- 
dure ll26l . The traversing parton gains a transverse momentum qx and radiates (be- 

fore or after the scattering) a gluon with a certain momentum k = (co, -^kx). The 
gluon differential distribution at first-order in opacity |60| is 

^JI^ ^ co^^- r^'d^q^ 2fcx-gx(fc-gx)^^" ^ ^28) 

where Xg is the mean free path of the radiated gluon. Applying the aforementioned 
recursive procedure, one obtains the gluon distribution to finite order (A^ > 1) in 
opacity. Each emission at a given opacity is assumed independent and a probabilistic 
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(b) (c) 
Fig. 9. Diagrams contributing to the lowest order in the opacity energy loss expansion fST]. 

scheme is set up wherein, identically to Eq. ( l26b , the parton loses an energy fraction 
e in n tries with a Poisson distribution 161], 



P„ie,E) 



-(m) 



-n"_ 



dl 
' d(i>i 



5(e£ 



7=1 



(29) 



where, (A^'^') is the mean number of gluons radiated per coherent interaction set. 
Summing over n gives the probability P(e) for an incident parton to lose a momen- 
tum fraction e due to its passage through the medium. This is then used to model a 
medium-modified FF, by shifting the energy fraction available to produce a hadron 
in a similar way as Eq. dZSl l. The key medium property to be obtained from the fits 
to the experimental data, is the initial gluon density dN^ /dy, after accounting for 
longitudinal expansion of the plasma. Note that the density of colour charges of a 
cylinder of plasma with "length" T and surface Ax, is p ~ dN^ /dy/{zA±). 

Higher Twist (HT) 

The higher-twist approximation ll83l [84l l85l l63l l64l l65l describes the multiple scat- 
tering of a parton as power corrections to the leading-twist cross section (Fig. [TOl ). 
These coiTections are enhanced by the medium-length L and suppressed by the power 
of the hard scale Q^. Originally, this approach was applied to calculate the medium 
connections to the total cross section in nuclear deep-inelastic eA scattering. 




Fig. 10. Higher-twist contributions to quark scattering in a medium (hatched area) 1861 . 



The scheme allows one to compute multiple Feynman diagrams such as those 
in Fig. [To] which are then combined coherently to calculate the modification of 
the fragmentation function directly as a medium-dependent additive contribution, 
with 



yarned 



'^vac 



fQ'dkia^ r^dx ^ r , 
Jo fc2 2% J,, X h.X 



J=<i,8 



)} 



(30) 
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Here, APi^j oc Pj^j Ca^X^T^^ represents the medium-modified splitting function of 
parton / into j (a momentum fraction x is left in parton j and the radiated gluon 
or quark carries away a transverse momentum k^^. The entire medium effects are 
incorporated in the nuclear quark-gluon correlation T^g term. The normalisation C of 
this correlator is set by fitting to one data point from which one can directly calculate 
the medium-modified FFs and then the final hadron spectrum. The parameter C can 
also be used to calculate the average energy loss suffered by the parton. 

AMY 

The Arnold-Moore-Yaffe (AMY) |l69l[87l[88llZl]l221 approach describes parton en- 
ergy loss in a hot equilibrated QGP, where a hierarchy T ^ gT ^ g^T is assumed. 
The hard parton scatters off other partons in the medium, leading to momentum trans- 
fers of ff{gT) and inducing collinear radiation. Multiple scatterings of the incoming 
(outgoing) parton and the radiated gluon are combined to get the leading-order gluon 
radiation rate. One essentially calculates the imaginary parts of ladder diagrams such 
as those shown in Fig.[TTl by means of integral equations which yield the 1^2 tran- 




sition rates F^^ of a hard parton a into a radiated gluon g and another parton b. These 
rates, with T-dependent Bose-Einstein (for gluons) and Fermi-Dirac (for quarks) 
exponential factors for the medium partons, are then used to evolve the original dis- 
tributions over the medium length by means of a Fokker-Planck like equation [71] 



dPcjp) 
dt 



(31) 



The medium modified FF are obtained from the convolution of the vacuum FF with 
the hard parton distributions when exiting the plasma \12\: 

= / dpf'^i^Pa{pf:pi)&:":,{7:). od 

J z „ 

where z — phj Pi and z' = ph/pf, with pi and pf the momenta of the hard partons 
immediately after the hard scattering and prior to exit from the medium respectively. 
The model of the medium is essentially contained in the space-time profile chosen 
for the initial temperature T appearing in the transition rates. 



Model comparison 

The four energy-loss formalisms discussed above can be roughly divided into two 
groups. Those calculating the radiated gluon spectrum i.e. the energy lost by the 
initial parton (BDMPS/ASW and GLV) and those determining directly the change 
in the final distribution of the traversing partons (Higher Twist and AMY). Each 
approach has its advantages and disadvantages: 
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• ASW and GLV: Both are applicable to thin and thick media (confined or decon- 
fined), but they do not account for the energy flow into the medium. 

• Higher Twist: It can directly compute the medium-modified fragmentation func- 
tions, and allows for the study of multi-hadron correlations, but the formalism is 
more appropriate for thin than thick media. 

• AMY: It is the only framework that accounts for processes where a thermal gluon 
or quark is absorbed by a hard parton, and elastic losses can be included in a sim- 
ple way, but it does not take into account vacuum radiation (nor vacuum-medium 
interference) and its application to non-thermalised media is questionable. 

All four schemes have independently made successful comparisons to the avail- 
able data (see Fig. fT2land coming Sections). The outcome of the models is one pa- 
rameter tuned to ideally fit all experimental observables: (q) in the BDMPS/ASW 
scheme, the initial dN^ /dy density in GLV, the C correlator normalization (or energy 
loss Eo) in HT, and the temperature T in AMY. All jet-quenching observables in AuAu 
collisions at ^^1^ = 200 GeV can only be reproduced with medium parameters con- 
sistent with a QGP at temperatures above the QCD phase transition (7^,,, « 0.2 GeV). 
The analytical results of the different schemes under "controlled" situations are in 
principle equivalent, see e.g. |7 |. Yet, the detailed comparison of the phenomenolog- 
ical results of the models is not always straightforward as they 

• use different approximations in their calculations (e.g. running-coupling, or fixed 
to a, = 0.3, 0.5), 

• do not always include the same list of physics processes (e.g. AEcn is neglected 
in some cases), 

• choose different fitting parameters to characterise the medium (see above), and 

• the space-time profile of the quenching medium is not always equivalent (e.g. 
static -plasma with average path-length vs. ID Bjorken expansion or vs. full 3D 
hydrodynamical evolutions with varying thermalisation times Xq). 
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Fig. 12. Suppression factor for high-pj- pions in central (top) and semicentral (bottom) AuAu 
collisions at RHIC I89J compared to AMY, HT and ASW energy loss calculations |90|. 



The quantitative consistency of the different schemes has been investigated 
within a 3-dimensional hydrodynamics approach (Fig. [T2T l ||90l Unking the various 
medium properties via thermodynamical relations and using the same space-time 
evolution. Yet, the extracted q values still differ by factors of 2 - 3 (see Section l4~2h ). 
At least part of the uncertainty is due to the relative insensitivity of the q parameter 
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to the irreducible presence of hadron from (unquenched) partons emitted from the 
surface of the plasma ||9T1 . Additional constraints on q can be placed by requiring 
also the model reproduction of the suppressed dihadron azimuthal correlations (see 
SectionlSTTTl. 

3.3 Jet quenching Monte Carlo's 

Ultimately, the energy loss schemes discussed in the previous Section are all based 
on a final energy-rescaling of the vacuum parton-to-hadron fragmentation func- 
tions (Fig. |7]). Recently, attempts to reformulate parton energy loss as a medium- 
modification of the perturbative evolution of the fragmentation functions have been 
implemented in Monte Carlo (MC) models [92 93 94 95 1. Such MC approaches 
allow one to address more detailed experimental observables such as the particle and 
energy flows within a jet. The DGLAP scale-dependence (Q^-evolution) equation of 
the FFs reads 



with splitting functions Pi^j{z), Eq. (fTOl i. The probabilistic nature of parton shower- 
ing- with a 5'M(iflA;ov/flcfor exp[— / dQ/Q^ J as/2KPi^j{z,Q^)dz] giving the prob- 
ability that a parton evolves from virtualities Qi to Q2 without branching - can be 
easily implemented in MC codes which compute the virtuaUty and energy fraction 
of a parton at each branching point with proper energy-momentum conservation. 
Parton showers are a basic ingredient of event generators such as PYTHIA |96| or 
HERWIG ||97ll often used to compare the experimental jet data to the details of the 
underlying QCD radiation pattern. Medium effects can be easily included by e.g. 
modifying the splitting functions in Eq. ( |33] |. HYDJET |98, 99J was the first MC 
code which incorporated medium effects via a PYQUEN routine which modifies the 
standard PYTHIA branching algorithm to include radiative and elastic energy losses. 
More recent developments like Q-PYTHIA and Q-HERWIG [MJ change the DGLAP 
evolution of these two parton-shower MCs. The JEWEL MC ||93l implements elastic 
scattering in DGLAP evolution plus radiative energy loss through a multiplicative 
constant in the infrared part of the splitting functions [33 1, whereas YaJEM [95) 
increases the virtuality of the traversing partons (i.e. their probabilities to radiate) in 
PYTHIA according to the medium q. 

3.4 Parton energy loss in AdS/CFT 

So far, we have discussed perturbative calculations of parton energy loss in an ideal 
QGP Yet, the medium produced at RHIC has temperatures ^{2T„i,) in a domain 
where lattice QCD [2J still predicts large deviations with respect to the asymp- 
totic ideal-gas behaviour. Many experimental signals at RHIC are consistent with 
the formation of a strongly-coupled plasma (sQGP) I1100l|10l1ll02| . Such a regime 
is theoretically treatable via the Anti-de-Sitter/Conformal-Field-Theory (AdS/CFT) 
correspondence between weakly-coupled gravity and strongly-coupled gauge theo- 
ries IfTOl . 

The AdS/CFT correspondence conjectures that string theories described in an 
Anti-de-Sitter space@ times a 5-dimensional sphere (AdSs x S5) are equivalent 
to a conformal field theory (CFT), defined on the 4-dimensional boundary of this 
space. A particularly useful case is ^ ~ 4 Super-Symmetric Yang Mills (SYmQ 

AdSs is a 5-dimensional space with constant and negative curvature. 
SYM is a quantum-field SU{Nc) theory like QCD = 4 indicates 4 additional super- 
charges) but dissimilar from QCD in many aspects: extra SUSY degrees of freedom, no 
running coupling, no confinement, ... Yet, such differences "wash out" at finite-J 11041 . 




(33) 
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at Strong coupling gym and large number of colours A^^. (i.e. at large 't Hooft coupling 
X = gym^c ^1) which is dual to supergravity in a curved space-time. The string cou- 
pling gs, the curvature radius R of the AdS metric, and the string tension (27ta')^' 
are related to the S YM quantities via: R^/a' = VX, and 4-%gs — gy^ — V-^c- Essen- 
tially, taking the large Nc limit at fixed A, (i.e. weakly coupled gravity: gs — > 0) and 
the large X limit (i.e. weakly curved space and large string tension), the SYM theory 
can be described by classical gravity in a 5-D space. By virtue of such a duality, 
one can carry out analytical calculations of gravity, which can then be mapped out 
"holographically" to the non-perturbative dynamics of the gauge (QCD-like) theory. 

One can further exploit the AdS/CFT correspondence for theories at finite- 
temperature, by replacing the AdSs space by an AdS Schwarzchild black-hole. The 
temperature of the gauge theory is then equal to the black-hole Hawking temper- 
ature: T = ro/inR'^), where ro is the coordinate of the black-hole horizon. Recent 
applications of this formalism in the context of heavy-ions physics have led to the 
determination of transport properties of strongly-coupled (SYM) plasmas - such 
as its viscosity 11_05|, the q parameter |106|, and the heavy-quark diffusion coeffi- 
cients II 1071 1 1081 II 091 ITTOl - from simpler black hole thermodynamics calculations. 

In the case of jet quenching calculations [iTT', '1 121, one can express the prop- 
agation of a parton through a medium in terms of Wilson lines. The q parameter is 
identified with the coefficient in the exponential of an adjoint Wilson loop averaged 
over the medium length: (W^(C)) °^ exp[^L] [106|. One then evaluates the grav- 
ity dual of this Wilson loop given by the classical action of a string stretching in 
an AdSs x Ss space with a Schwarzschild black hole background. After solving the 
equations of motion of the string, the q parameter is found to be 



Though this result is computed in the infinite coupling and number of colours limits, 
typical values of a, = 0.5 and A^,. = 3 lead to ^ = 4.5 - 20.7 GeV^/fm for T = 0.3 - 
0.5 GeV 111 121 . consistent with phenomenological fits of the RHIC data ||9T1 . There 
have been also AdS/CFT-based calculations 1 107] \lM HHH [TTOl of the transport 
properties of a heavy quark, described by a semiclassical string in the gravity theory, 
such as its diffusion constant in a ^ = 4 SYM plasma II108I 



which agrees with the drag coefficient, see Eq. (O, computed independently. 
4 High- pt leading hadron suppression: data vs. theory 

The simplest empirically testable (and theoretically computable) consequence of jet 
quenching in heavy-ion collisions is the suppression of the single inclusive high- 
Pt hadron spectrum relative to that in proton-proton collisions. Since most of the 
energy of the fragmenting parton goes into a single leading hadron, QCD energy loss 
was predicted to result in a significantly suppressed production of high-pj- hadrons 
(Raa ^1) EH- We compare in this Section the existing measurements of large- 
Pt hadroproduction in pp and AA collisions, and discuss their agreement with jet 
quenching models. 



^sym — 




(34) 




(35) 
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4.1 High- pt hadron spectra in proton-proton and proton-nucleus collisions 

Figure [T3]collects several pj- -differential inclusive cross sections measured at RHIC 
in pp collisions at y/s = 200 GeV: jets Ml 13L charged hadrons Ml 14L neutral pi- 
ons Ml 151 . direct photons [ 116|, and D,B mesons (indirectly measured via inclusive 
from their semileptonic decays) |117| at central rapidities (y — 0), and nega- 
tive hadrons at forward pseudorapidities (r| = 3.2) Bl 181 . The existing measurements 
cover 9 orders of magnitude in cross section (from 10 mb/GeV^ down to 1 pb/GeV^), 
and broad ranges in transverse momentum (from zero for D,B mesons up to 45 
GeV/c, a half of the kinematical limit, for jets) and rapidity (rj = - 3.2). 




Fig. 13. Compilation of hard cross sections in pp at ^/s = 200 GeV measured by STAR 11131 
[TT4l . PHENIX ||TT5][Tl6][TlIl, and BRAHMS GH) (10%-30% syst. uncertainties not shown 
for clarity) compared to NLO I119II1201 and NLL 11211 pQCD predictions (yellow bands). 

Standard next-to-leading-order (NLO) 111 191 11201 or resummed next-to-leading 
log (NLL) |ii2r| pQCD calculations (yellow bands in Fig. \T3[ with recent proton 
PDFs II122L fragmentation functions 11231 1124L and with varying factorisation- 
renormalisation scales (/j = pr /2 — Ipr) reproduce well the pp data. This is true 
even in the semi-hard range pj ^ 1—4 GeV/c, where a perturbative description 
would be expected to give a poorer description of the spectra. These results indi- 
cate that the hard QCD cross sections at RHIC energies are well under control both 
experimentally and theoretically in their full kinematic domain. 

Not only the proton-proton hard cross sections are well under theoretical control 
at RHIC but the hard yields measured in deuteron-gold collisions do not show either 
any significant deviation from the perturbative expectations. Figure [14] shows the 
nuclear modification factors measured in dAu collisions at y^s^ = 200 GeV for high- 
er n° at y = iQlSllTia. The maximum deviation from the RdAu = 1 expectation 
is of the order of ~10%, well accounted for by standard pQCD calculations IJ271 
11281 that include DGLAP-based parametrisations of nuclear PDFs II 1291 and/or a 
mild amount of initial-state pj broadening [130] to account for a modest "Cronin 
enhancement'El 11311 . [The only exception to this behaviour is baryon (in particular. 

The "Cronin effect" is the observation of enhanced hadron production at ~ 1 -7 GeV/c 
in proton-nucleus compared to proton-proton collisions (see L130J for a recent review). 
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Fig. 14. Nuclear modification factors for high-pr Jt" (left) and r| (right) mesons at midrapidity 
in dAu collisions at ^/i^ = 200 GeV I125II126I . 



proton) production, which shows a large Cronin enhancement: Riiau = 1.5-2.0 f 1321.] 
These data clearly confirm that around midrapidity at RHIC energies, the parton flux 
of the incident gold nucleus can be basically obtained by geometric superposition 
of the nucleon PDFs, and that the nuclear {x,Q^) modifications of the PDFs are 
very modesQ Since no dense and hot system is expected to be produced in dAu 
collisions, such results indicate that any value of R^a different than sa 1 ± 0.2 
potentially observed for hard probes in AuAu collisions (at central rapidities) can 
only be due to final-state effects in the medium produced in the latter reactions. 

4.2 High-pT hadron spectra in nucleus-nucleus collisions 

Among the most exciting results from RHIC is the large high-pj- hadron suppression 
{Raa 1) observed in central AuAu compared to pp or dAu reactions. We discuss 
here the properties of the measured suppression factor and compare it to detailed 
predictions of parton energy loss models. 

(a) Magnitude of the suppression: medium properties 

Figure [15] shows the T\P spectrum measured in pp collisions II 1341 compared to pe- 
ripheral (left) and central (right) AuAu spectra f89l| at 200 GeV, as well as to NLO 
pQCD calculations 111 19L Whereas the peripheral AmAm spectrum is consistent with 
a simple superposition of individual A^A^ collisions, the data in central AmAm show a 
suppression factor of 4 - 5 with respect to this expectation. The amount of suppres- 
sion is better quantified taking the ratio of both spectra in the nuclear modification 
factor, Eq. (|23] |. Figure[T6]compiles the measured Raa [pt) for various hadron species 
and for direct y in central AuAu collisions at ^fs^ = 200 GeV. 

Above pT^5 GeV/c, 71° [1191, r| |135|, and charged hadrons lfTT4l[T40l (dom- 
inated by Jt* II140I ) show all a common factor of ~5 suppression relative to the 
Raa = 1 expectation that holds for hard probes, such as direct photons, which do no 
interact with the medium 1I136L The fact that Raa ~ 0.2 irrespective of the nature 
of the finally produced hadron is consistent with a scenario where final-state energy 
loss of the parent parton takes place prior to its fragmentation into hadrons in the 
vacuum according to universal (but energy-rescaled) FFs. The suppression factor at 

The same is not true a.i forward rapidities where gluon saturation effects in the Au PDFs 
play an important role (see e.g. 11331 and references therein). 
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Fig. 15. Invariant 7t" yields measured by PHENIX in peripheral (left) and central (right) 
AuAu collisions (squares) [89|, compared to the (J^A-scaled) pp ^ -\- X cross section (cir- 
cles) 11341 and to a NLO pQCD calculation (curves and yellow band) 11191 . 
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Fig. 16. Raa{pt) measured in central AuAu at 200 GeV for Jl" t89l and T| 11351 mesons, 
charged hadrons f 1 141, and direct photons I136II1371 compared to theoretical predictions for 
parton energy loss in a dense medium with dN^ /dy = 1400 (yellow curve) 11381 . 



top RHIC energies is very close to the "participant scaling", {Npan /2) /Ncoii ~ 0. 17, 
expected in the strong quenching limit where only hadrons coming from partons 
produced at the surface of the medium show no final-state modifications in their 
spectra II141I . From the Raa one can approximately obtain the fraction of energy 
lost, e/oj., = I^ptIpt, via 

e/„.v.v«l--R^*""^\ (36) 

when the AuAu and pp invariant spectra are both a power-law with exponent «, i.e. 
1 /prdN/dpT °^ pj" El. At RHIC (n « 8, /?a4 « 0.2), one finds £/„,,. w 0.2. 

The high- pt AuAu suppression can be well reproduced by parton energy loss 
models that assume the formation of a very dense system with initial gluon ra- 
pidity densities dN^ /dy « 1400 (yellow line in Fig. [TSI l II138L transport coeffi- 
cients {q) fa 13 GeV^/fm (red line in Fig. [17] left) |78|, or plasma temperatures 
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T w 0.4 GeV [72j. The quality of agreement between the theory and data has been 
studied in detail in ||89l 11431 taking into account the experimental (though not theo- 
retical) uncertainties. The PHENIX Ti" suppression data allows one to constrain the 
transport coefficient of the PQM model |78| {q) as 13.2 +32 and +52 GeV^/fm at the 
one and two standard-deviation levels (Fig. [17] right). 
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Fig. 17. Left: Raa{pt) for neutral pions in central AiiAu collisions (triangles) f89l com- 
pared to PQM predictions |78 | for varying values of the q coefficient (red curve, best fit 
for (q) = 13.2GeV^/fm). Right: Corresponding (data vs. theory) values for the PQM q 
parameters that fit the data points on the left plot H43j . 



The consistency between the extracted q, dN^ /dy and T values in the various 
models can be cross-checked considering the simple case of a gluon traversing a 
thermalised gluon plasma. The transport coefficient, Eq. (|2|i, is the product of the 
medium particle density, the medium Debye-mass, and the parton-medium cross 
section. Taking <3gg = 97ta^/(2m|,) with a, = 0.5 for the latter, one has a simple 
relatior0 between q and p: 

^[GeV2/fm] = X O X p = x 97ia^/(2OT|,) x p « Q.UKp[fm-^] , (37) 

where we introduce a /T-factor to account for possible higher-order scatterings not 
included in the LO perturbative expression for Ogg. For an ideal ultrarelativistic gas, 
the particle density scales with the cube of the temperature as p w ndf/9 • T^. For a 
pure gluon plasma (with ndf = 16 degrees of freedom), p[fm"^] w 260 • (r[GeV])^, 
and one can write Eq. ( l37T i as: 

^[GeV^/fm] « 36/: • (7'[GeV])^ (38) 

In addition, from the relation p[fm"^] w 1.9 • (e[GeV/fm^])^/^ between particle and 
energy densities, one can also express Eq. (|37| i as: 

^[GeV^/fm] w 0.21 K- (eLGeV/fm^])^/^ . (39) 

In an expanding plasma, the density follows a power-law evolution as a function of 
time, p = po ('Co/'c)", and thus so does the transport coefficient (|37] |: 

«GeV2/fm]«0.14ir.pom" = 0.14/r.^n«0.14/r.-l^l, 

V T / dV \ z J At dy 1 

(40) 

Conversion between units is done multiplying by suitable powers of he = 0. 197 GeV fm. 
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where for the two last equahties we have assumed a 1 -dimensional (aka Bjorken) 
longitudinal expansion i.e. a = 1 and dV = At Iq dy, where A7'[fm^] is the transverse 
area of the system. Combining Eq. ( |40] i with Eq. jZTl ) that relates the time-averaged 
q{x) to that of a static medium with effective length L^ff, we finally get 

2 dN^ o dN^ 

(q) [GeV2/fm] « 0. 14^: • . « 1 .4 IQ-^ • • — , (41) 

Leff[tm]AT[fnr] dy dy 

where, for the last equality, we use L^jf ~ 2 fm and {At) ~ 100 fm^ for the overlap 
area in 0-10% most central AuAu. This approximate relation between the average 
transport coefficient and the original gluon density is only well fulfilled by the data 
(see Table 12] below) for very large K 1 factors. The fact that the jet-quenching 
data favours an effective elastic parton-medium cross-section much larger than the 
LO perturbative estimate of <3gg 1.5 mb, see Eq. dJTl l. has been discussed many 
times in the literature - e.g. in the context of the strong partonic elliptic flow seen 
in the data Ml 441 - and supports the strongly-coupled nature of the QGP produced at 
RHIC |fT02l. 



Equation is just a simple order-of-magnitude estimate based on simplifying 
assumptions. A more detailed comparison of different energy-loss schemes within a 
realistic 3 -dimensional hydrodynamics evolution has been carried out in [145J . The 
extraction of a common q parameter from the different model predictions relies on 
the use of (thermo)dynamical relationships such as Eqs. ( |38] | or ( |39] l. The results for 
the ASW, AMY and HT schemes are shown in Table [T] The ASW calculations con- 
sistently predict a higher q than AMY or HT. Seemingly, as of today, comparisons of 
model predictions to RHIC results for Raa{pt) can only constrain q within a factor 
of 2 - 3. The origin of such a large variability can be traced to a combination of 
(i) the relative insensitivity of using just a single-inclusive observabl^^ Raa{pt), 
in the data-model comparisons 191] (additional independent measurements place 
extra constraints on q as discussed in Section 15.1b . and (ii) the assumptions about 
the equation-of-state of the medium (and its time evolution) and the correspond- 
ing approximations relating its thermodynamical and transport properties. Genuine 
model differences (e.g. AMY accounts for collisional loses which are neglected in 
the purely radiative ASW approach) play also a role. A working group 01461 has 
been recently created to clarify discrepancies among the formalisms. 



Table 1. Transport coefficients q derived in a 3-D hydro simulation of an expanding QGP with 
initial temperature Tq = 0.4 GeV (at To = 0.6 fm/c) 11451 with different parton energy loss 
implementations (ASW, HT and AMY schemes) that reproduce the high-pj- suppression 
observed in central AmAm at RHIC [89 1. The a,h exponents indicate two choices of scaling of 
^(r,T) with the initial plasma temperature or energy-density: (a) °= ^^'(•'i'^)' ''nd {h) go °= 
e0'''*(r, t). The PQM/ASW result (Fig. [17] (q) for a static plasma) is also listed for comparison. 



ASW 



HT 



AMY 



q (GeV^/fm) lO'"' - IS.sC'), 13.2(^2^) 2.3('') -4.3(*) 4.l(") 



Irreducible parton production from the outer corona of the medium - which remains un- 
suppressed even for extreme densities in the centre - makes of Rjy^ (pT ) a "fragile" observ- 
able [91J. 
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(b) Centre-of-mass energy dependence 

As one increases the centre-of-mass energy in nucleus-nucleus collisions, the pro- 
duced plasma reaches higher energy and particle densities and the system stays 
longer in the QGP phase. Since AEi^ss °= dN^ /dy oc dNch/dr[, and since the charged 
particle multiplicity in AA at midrapidity increases with collision-energy as [147 J 

dN,h/dy\ « 0.75 • {N,a„/2) ■ ln( ^AGeV]/1.5), (42) 

(where Npan is the number of nucleons participating in the collision), one naturally 
expects the hadron quenching (at a given centrality) to increase logarithmically with 
y^J^. The actual "excitation function" of Raa does not follow exactly the same c.m.- 
energy dependence of Eq. (|42j because for increasing energies other factors play 
counteracting roles: (i) the lifetime of the quenching medium becomes longer (which 
enhances the energy loss), (ii) the parton spectrum becomes flatter (which leads to 
a comparatively smaller suppression for the same value of AEioss, see below), and 
(iii) the relative fraction of quarks and gluons produced at a given pj changes and so 
does the quenching factor (see Fig.|22]below and the discussion on the colour-factor 
dependence of the suppression). 



High tt" in 0-10% central A+A: 
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AuAu \fs;;;;;= 62.4 GeV (PHENIX prelim.) 
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Fig. 18. Nuclear modification factor, Raa{Pt), for neutral pions in central PbPb at 
= 17.3 GeV 1 141, i49J and AuAu at = 62.4 GeV QSOJ, 200 GeV [139]; com- 

pared to GLV energy loss calculations for initial gluon densities: dN^ /dy = 400, 800, 1400 
I138II1511 respectively. Experimental normalisation errors, i^7'(10%-25%), not shown. 



Figure [TS] compiles the measured Raa{pt) for high-y^j- 71'' measured in central 
AA collisions in the range ^/s^ w 20 - 200 GeV compared to parton energy loss 
calculations that assume the formation of a QGP with initial gluon densities per unit 
rapidity in the range dN^ /dy w 400 - 1400 [138, 151] or, equivalently, averaged 
transport coefficients (q) w 3.5 - 13 GeV^/fm llTSl (Table |2]l. The tiieoretical pre- 
dictions reproduce very well the experimental data. The SPS data show an Raa for 
central PbPb which, though consistent with unity [149|, is significantly suppressed 
compared to the "Cronin enhancement" observed for peripheral PbPb and for pPb 
collisions fT52l. The onset of the suppression must lie close to the highest ener- 
gies reached at SPS, in a domain of cm. energies that would be worth studying 
detailedly (including proper high-/? 7- pp reference measurements) in coming RHIC 
runs. At y^s^ = 62 GeV, the suppression is already comparatively large because, as 
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can be seen from Eq. ( [36] ), Raa not only depends on A£/„„ but also on the steep- 
ness (power-law exponent n) of the parton pj spectrum: with decreasing ,Js^, the 
Pt spectra become steeper effectively leading to a comparatively larger suppression 
(i.e. smaller Raa) for the same value of AEioss- 

In any case, it is interesting to remark that for each collision energy the val- 
ues for dN^ /dy derived from the jet-quenching models are consistent with the final 
charged hadron density dNch/dr\ measured in the reactions. This is expected in an 
isentropicF^ expansion process, where all the hadrons produced at midrapidity in a 
AA collision come directly from the original gluons releasee^ in the collision; 



dy ^ Nch 



dy\ 

dy 



^^1.8.^. (43) 
dr\ d^ 



This relationship is relatively well fulfilled by the data as can be seen by comparing 
the fourth and fifth columns of Table |2] 



Table 2. Initial gluon densities dN^ /dy II138II1511 , and transport coefficients {q) r78] |153l 
for the dense media produced in central AA collisions at SPS and RHIC obtained from parton 
energy loss calculations reproducing the observed high-pj- Jl^ suppression at each ^/s^. The 
measured charged particle densities at midrapidity, dN^ /df] L147J . are also quoted. 





(GeV) 


(GeV^/fm) 


dNS/dy 




SPS 


17.3 


3.5 


400 


312 ±21 


RHIC 


62.4 


7. 


800 


475 ± 33 


RHIC 


130. 


~11 


~1000 


602 ± 28 


RHIC 


200. 


13 


1400 


687 ± 37 



(c) -dependence of the suppression 

At RHIC top energies, the hadron quenching factor remains relatively constant from 
5 GeV/c up to the highest transverse momenta measured so far, pT ~ 20 GeV/c (see 
Figs.[T6l[T9l). On rather general grounds 1 154|, one expects a rise of Raa with pj for 
any model in which the energy loss probability does not strongly depend on the initial 
parton energy as more of the shift in energy becomes accessible. The detailed form of 
the rise is sensitive to the energy loss probability distribution, Eq. (l26T l. The measured 
flatness of Raa{pt) was not expected in various original analytical QCD energy-loss 
calculations including the LPM effect (see e.g. II 156 1) which instead predicted an 
Raa slowly (logarithmically) increasing with pr- However, the combined effect of 
(i) kinematics constraints (which modify the asymptotic AEi^ss formulas), (ii) the 
steeply falling pr spectrum of the scattered partons, and (iii) i^(20%) pj- -dependent 
(anti)shadowing differences between the proton and nuclear PDFs included in the 
various models [138. J8.i71ii91 j , do result in an effectively flat RAAipr) as found in 
the data. 

The much larger kinematical range opened at LHC energies [155] will allow one 
to test the pj- -dependence of parton energy loss over a much wider domain than at 
RHIC. The GLV and PQM predictions for the charged hadron suppression in PbPb 
at 5.5 TeV are shown in Fig. [19] Apart from differences in the absolute quench- 
ing factor, PQM seemingly predicts a slower rise of Raa{pt) than GLV. The large 

Namely, expanding at constant entropy i.e. without extra particle production. 

We use: Niot /Nch = 3/2 and the Jacobian Idti/rfyl =E/mj ~ 1.2 for a mostly pionic system. 
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Fig. 19. Raa{pt) for neutral pions at SPS 1148111491 and RHIC li89J compared to the suppres- 
sion of charged hadrons in central PbPb at the LHC ( = 5.5 TeV) predicted by the GLV 
{dm /dy = 2000 - 4000) (l38l[T55] and PQM {{q) 30 - 80 GeV^/fm) (HIlSS] models. 



PT reaches of the ALICE IT5T\, ATLAS lfT58l and CMS |fT59l experiments (up to 
300 GeV/c for the nominal luminosities) will allow them to test such level of model 
details. 



(d) Centrality (system-size) dependence 

The volume of the overlap zone in a heavy-ion collision can be "dialed" either by 
selecting a given impact-parameter b - i.e. by choosing more central or peripheral 
reactions - or by colliding larger or smaller nuclei. From Eq. (fT4l i. the relative amount 
of suppression dependo on the effective mass number A^s or, equivalently, on the 
number of participant nucleons in the collision A^pait, as: e — AE/E 
Combining this expression with Eq. ( l36l l yields II 1601 

= (1 - KA^"„rt)""^ J with a w 2/3, and k an arbitrary constant. (44) 

Figurel20l(left) compares the measured high-pj- pion suppression in CuCu and AmAm 
at y/s^ = 200 GeV II161II162L Because of the large difference in the Cu (A = 63) 
and Am (A = 197) mass numbers, the sameNpart value (i.e. the same overlap volume) 
implies very different collision geometries: a thin, elongated collision zone in AuAu, 
and a thicker, more spherical one in the CuCu case. Yet, interestingly the average 
suppression in the two systems depends only on Npari- Fitting this dependence with 
expression ( l44b yields a — 0.56 ± 0. 10, consistent with a w 0.6 expected in detailed 
parton energy loss calculations fW TSOl. 

The right plot of figure |20l shows the Raa{pt) measured in CuCu at 22.4, 62.4, 
and 200 GeV II161I . The observed amount of suppression is roughly a factor of 
(Aam/Ach)^/^ w 2 lower than in AuAu at the same energies (Fig. [TSll. The Raa{pt) 
can be described by the GLV model with initial gluon densities dN^ / dy^ 100 - 370 
(the CuCu enhancement at 22.4 GeV is actually consistent with a scenario without 
parton energy loss). 



(e) Path-length dependence 

The analytical quadratic dependence of the energy loss on the thickness of a static 
medium L, Eq. ( fT2] l. becomes effectively a linear dependence on the initial value of 

3' Sincedm/dy oc dN^h/dy oc A^ff oc A^p^^, LocA^^ oc N^i^^, andA^ocA%^ oc N^l^. 
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Fig. 20. Left: Centrality (Npart) dependence of the high-pj- Jl" suppression in CuCii and AkAm 
at 200 GeV 1 1 62 1 . Riglit: Rm {pt ) for 7l" in central CuCu collisions at 22.4, 62.4 and 200 GeV 
compared to GLV calculations with initial gluon densities dN^ /dy ~ 100 - 370 [ 161 j . 



L when one takes into account the expansion of the plasma, see Eq. ( fT4l ). Experimen- 
tally, one can test the L-dependence of parton suppression by exploiting the spatial 
asymmetry of the system produced in non-central nuclear collisions (Fig. [21] left). 
Partons produced "in plane" ("out-of-plane") i.e. along the short (long) direction of 
the ellipsoid matter with eccentricity e will comparatively traverse a shorter (longer) 
thickness. 




Fig. 21. Left: Effective thicknesses along various azimuthal directions with respect to the 
reaction plane in the overlap region of two heavy-ions. Right: Fraction of energy loss e/oj^ 
versus effective path-length measured for high-pj neutral pions in AuAu at 200 GeV 11421 . 



PHENIX 01421 11631 has measured the high-/?7- neutral pion suppression as a 
function of the angle with respect to the reaction plane, Raa{pt,<^)- Each azimuthal 
angle (|) can be associated with an average medium path-length Lg via a Glauber 
model. Figure|2T](right) shows the measured fractional energy loss e/oii((|)), obtained 
via Eq. ( |36] |. as a function of Lg for pions in the range pj- = 5 - 8 GeV/c (markers of 
different colours coiTespond to varying centralities, i.e. eccentricities e). The energy 
loss is found to satisfy the expected AEu,ss °= L dependence above a minimum length 
of L ~ 2 fm. The absence of suppression in the surface of the medium is explained 
as due to a geometric "corona" effect 1 164|. 
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(f) Non-Abelian (colour factor) dependence 

The amount of energy lost by a parton in a medium is proportional to its colour 
Casimir factor Cr, i.e. Ca = 3 for gluons, Cp = 4/3 for quarks. Asymptotically, the 
probability for a gluon to radiate another gluon is Ca/Cf = 91 A times larger than for 
a quark and, thus, g-jets are expected to be more quenched than ^-jets in a QGP. One 
can test such a genuine non-Abelian property of QCD energy loss in two ways: 

(1) by measuring hadron suppression at a. fixed pj for increasing ^/s II153II1651 . 

(2) by comparing the suppression of high-pj- (anti)protons (coming mostly from 
gluon fragmentation) to that of pions (which come from both g and q^q). 

The motivation for (1) is based on the fact that the fraction of quarks and gluons 
scattered at midrapidity in a pp or AA collision at a fixed pr varies with y^s^ in a 
proportion giveioby the relative density of q,q and g at the corresponding Bjorken 
X — 2pT / "v/s in the proton/nucleus. At large (small) x, the hadronic PDFs are dom- 
inated by valence-quarks (by "wee" gluons) and consequently hadroproduction is 
dominated by quark (gluon) scatterings. A full NLO calculation [1191 (Fig. [221 left) 
predicts that hadrons with pr ^ 5 GeV/c at SPS (LHC) energies are ~100% pro- 
duced by quarks (gluons), whereas at RHIC they come 50%-50% from both species. 




Fig. 22. Left: Relative fraction of quarks and gluons fragmenting into a hadron at pj- = 5 GeV/c 
in pp collisions in the range ^ = 10 - 5500 GeV given by NLO pQCD H19J . Right: 
Raa{Pt = 4 GeV/c) for ifi in central AA collisions as function of collision energy compared 
to non-Abelian (solid) and "non-QCD" (dotted) energy loss curves I153II1651 . 



Figure |22] (right) shows the Raa for 4-GeV/c pions measured at SPS and RHIC 
compared to two parton energy loss curves, both normalised at the Raa ~ 1 measured 
at SPS and extrapolated all the way up to LHC energies [1651. The lower curve shows 
the expected Raa assuming the normal non-Abelian behaviour (AEg/AEq = 9/4). The 
upper (dotted) curve shows an arbitrary prescription in which quarks and gluons lose 
the same energy (AEg = AEg). Above ^/s^ « 100 GeV, gluons take over as the dom- 
inant parent parton of hadrons with pr ^ 5 GeV/c and, consequently, the Raa values 
drop faster in the canonical non-Abelian scenario. The experimental high-pj- Tt" data 
thus supports the expected colour-factor dependence of Raa ( ^/s^) [1531 ■ 

The second test of the colour charge dependence of hadron suppression is based 
on the fact that gluons fragment comparatively more into (anti)protons than quarks 
do. One would thus naively expect R'^ < R^^. The STAR results (Fig. |23l left) 



The different "hardness" of quarks and gluons fragmenting into a given hadron at the cor- 
responding z = Phadron/P parton plays also a (smaller) role. 
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are however seemingly at variance with this expectation: pions appear more sup- 
pressed than protons at high-pj- |I166|. The use of (anti)protons as a reference for 
perturbative particle production is however questionable: p,p are already found to 
be enhanced in dAu compared to pp collisions by a factor ~ 50% - 100% for pr's 
as large as 7 GeV/c 1 132|. It is likely that there is an extra mechanism of baryon pro- 
duction, based e.g. on in-medium quark coalescence 11681 . which compensates for 
the energy loss suffered by the parent partons. It is also important to stress that the 
AEg /AEg = 9/4 expectation holds only for asymptotic parton energies. Finite energy 
constraints yield AEg/AEq w 1.5 for reaUstic kinematics (Fig.|23] right) l,78.,167J . 

(g) Heavy-quark mass dependence 

A robust prediction of QCD energy loss models is the hierarchy AEq < AEq < AEg. 
Due to the dead-cone effect, the radiative energy loss for a charm (bottom) quark is 
--25% (75%) less than for a light-quai-k (see Section Ell- Surprisingly, PHENIX 
and STAR measurements of high-/?/' electrons from the semi-leptonic decays of 
D- and B-mesons (Fig. |24] | indicate that their suppression is comparable to that of 
light mesons: Raa{Q) ^ RAA{q,g) ~ 0.2 |169, 170l, [TtTI . Such a low Raa cannot 
be described by radiative energy loss calculations with the same initial gluon den- 
sities or transport coefficients needed to reproduce the quenched light hadron spec- 
traimiimi. 

Various explanations have been proposed to solve such a 'heavy flavor puzzle' : 

• First, if only c quarks (three times more suppressed than the heavier b quarks) 
actually contributed to the measured high-pj- decay electron spectrum, then 
one would indeed expect ^a4(c) « 0.2 ITTTI . Yet, indirect measurements from 
PHENIX 1 175 1 and STAR [1761 have confirmed the similar production yields of 

from D and B decays above pj ^ 5 GeV/c predicted by NLL pQCD |.177j . 

• The heavy-quark suppression has revived the interest of computing elastic energy 
loss in a QGP 1 17l[n9l|l9l|T80l|20l. As discussed in Sect.O AE.^ii can indeed 
be a significant contribution for heavy quarks (see 'rad.H-coll.' curves in Fig.l24l 
left). 

• The strongly-coupled nature of the plasma at RHIC would lead, according to 
AdS/CFT calculations jW] [Ml \M\ \lM [Ml, to a larger heavy-quark mo- 
mentum diffusion parameter than expected in perturbative approaches II 1821 . This 
could explain the larger charm/bottom quenching observed in the data. 

• Two works [183. ,184J have argued that the large charm-quark coalescence into 
Ac baryons (with a small semileptonic decay branching ratio) in the plasma 
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Fig. 24. Raa{pt) for decay electrons from D and B mesons in central AuAu at y'SiviV = 
200 GeV |169l|170lll71| compared to various radiative+elastic energy loss models for c and 
b quarks (left) and to a model of D and B meson dissociation in the plasma 11741 (right). 



would deplete the number of open-charm mesons, and correspondingly reduce 
the number of decay electrons, compared to pp collisions. 
The assumption of vacuum hadronisation (after in-medium radiation) implicit in 
all parton energy loss formalisms may well not hold in the case of a heavy quark. 
All existing quark-hadronisation time estimates B 13011 are inversely proportional 
to the mass m/, of the final produced hadron: the heavier the hadron, the fastest it 
is formed. In the rest frame|2j of the fragmenting heavy-Q, the formation time of 
D- and B-mesons II 1741 

V™. = tTr- ^.2^n ^'^^^I'K — ' ^^^'^ = pq/Eq, (45) 

is of order T^^^.^^_ « 0.4 - 1 fm/c respectively. Thus, theoretically, one needs to 
account for both the energy loss of the heavy-quark as well as the possible dis- 
sociation of the heavy-quark meson inside the QGP. The expected amount of 
suppression in that case is larger and consistent with the data (Fig.|24l right). 



5 High-pr di-hadron (|),ri correlations: data vs. theory 

Beyond the leading hadron spectra discussed in the previous Section, detailed studies 
of the modifications of the jet structure in heavy-ion collisions have been addressed 
via high- Pt multi -particle (mostly di-hadron) (|),ri correlations at RHIC (and to a 
lesser extent at SPS). Jet-like correlations are measured on a statistical basis by se- 
lecting high-pj- trigger particles and measuring the azimuthal (A(|) = (|) — i^n-ig) and 
pseudorapidity (Ar| = r| — ri,,-,^,) distributions of associated hadrons (p^""' < pj'^) 
relative to the trigger: 

C(A(^,Aii) = — (46) 

Combinatorial background contributions, corrections for finite pair acceptance, and 
the superimposed effects of collective azimuthal modulations (elliptic flow) can 
be taken care of with different techniques ||185[ 11861 11871 . A commonly-used 
C(A(|)) background-subtraction method is the so-called "zero yield at minimum" 
(ZYAM) lfT88l . 



Note that in the laboratory system there is an extra Lorentz boost factor: x,,^^ = Jq ■ 
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Fig. 25. Schematic illustration summarising the jet-induced di-hadron correlation signals in 
A(|) and At] observed in pp (left) and central AuAu (right) at = 200 GeV L189J . 

A schematic representation of the di-hadron azimuthal-pseudorapidity coiTe- 
lations dNi,air/<iA(^dAr[ measured in pp and central AuAu collisions is shown in 
Fig.|25] In the pp case, without significant initial- or final- state interactions, a dijet 
signal appears clearly as two distinct back-to-back Gaussian-like peaks at A(|) « 0, 
Ar| w (near-side) and at A(|) w 7i (away-side). Note that the away-side peak is nat- 
urally broader in Ar| (up to Ar| w 2) than the near-side peak due to the longitudinal 
momentum imbalance between the two colliding partons with different xi,X2 mo- 
mentum fractions (the collision is boosted in r| by an amount ln(xi /x2) in the direc- 
tion of the larger jici,2). At variance with such a standard dijet topology, the di-hadron 
correlations in AuAu reactions at RHIC show several striking features, discussed in 
detail below: 

• The away-side azimuthal peak at A(|) « 7t is strongly suppressed with increasing 
centrality for hadrons with pf "'^' > 2 GeV/c, consistent with strong suppression 
of the leading fragments of the recoiling jet traversing the medium 111 851 . 

• The vanishing of the away-side peak is accompanied with an enhanced produc- 
tion of lower pt hadrons {pj™'^' < 2 GeV/c) LI 86 J with a characteristic double- 
peak structure at A(|) w 7i± 1.1 - 1.3 lfT86l[T87l . 

• One observes a large broadening ("ridge"), out to Ar| « 4, of the near-side 
pseudo-rapidity c/A^pair/t/Ari correlations l,187J . 



5.1 Azimuthal correlations: away-side quenching and energy loss 

Figure |26] shows the increasingly distorted back-to-back azimuthal coiTelations in 
high-/?7- triggered central AuAu events as one decreases the pr of the associated 
hadrons (right to left). Whereas, the AuAu and pp near-side peaks are similar for all 
Pt's, the away-side peak is only present for the highest partner pr's but progressively 
disappears for less energetic partners II190II1911 . Early STAR results 01851 showed a 
monojet-like topology with a complete disappearance of the opposite-side peak for 
possoc ~ 2 - 4 GeV/c. 

For any range of trigger p'-f'^ and associated p'^""^' intervals, the correlation 
strength over an azimuthal range A(|) between a trigger hadron /i, and a partner hadron 
ha in the opposite azimuthal direction can be constructed as a function of the momen- 
tum fraction zj- — pj'"'"^'/pj'^ via a "pseudo-fragmentation function" II192I : 



^AA y^T)- dp J dp J / dA<^ 



(47) 
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Fig. 26. Comparison of the azimuthal di-hadron correlation dNpair /dA(^dr[ for pp (open sym- 
bols) and central AuAu (histograms) at = 200 GeV for p'"^ = 5-10 GeV/c and increas- 
ingly smaller (right to left) values of pf""^ 11911 . 



Figure|22](left) shows the measured D"J^^ distributions for p p and AuAu collisions as 
a function of zt compared to predictions of the HT jet-quenching model for various 
values of the Eq parameter quantifying the amount of parton energy loss 1,193 1 . Sim- 
ilarly to Raa, the magnitude of the suppression of back-to-back jet-like two-particle 
correlations can be quantified with the ratio Iaa{zt) = Daa{zt)/Dpp{zt)- I^a^ (^t) 
(bottom-left panel of Fig.lZTli is found to decrease with increasing centrality, down to 
about 0.2 - 0.3 for the most central events II185II1941 . The right plot of Fig.|27]shows 
the best Eq w 1.9 GeV/fm value that fits the measured Raa and Iaa factors. Due to the 
ureducible presence of (unquenched) partons emitted from the surface of the plasma, 
the leading-hadron quenching f&ctOT Raa{pt) is in general less sensitive to the value 
of Eo than the dihadron modification ratio Iaa{zt)- The combination of Raa{pt) and 
Iaa{zt) provides robust quantitative information on the medium properties. 




Fig. 27. Left: 0"^^-" {zt) distributions for dAu and AuAu collisions at 200 GeV and Iaa{zt) 
ratio (for central AuAu) 1 194], compared to HT calculations 1 193 1 for varying Eq energy loss 
parameter. Right: Corresponding (data vs. theory) values for the fitted Eq parameters 11931 . 



5.2 Azimuthal correlations: away-side broadening and "conical" emission 



Since energy and momentum are conserved, the "missing" fragments of the away- 
side (quenched) parton at intermediate pj-'s should be either shifted to lower energy 
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(Pt ^ 2 GeV/c) and/or scattered into a broadened angular distribution. Both, soft- 
ening and broadening, are seen in the data when the pj of the away-side associated 
hadrons is lowered (see two leftmost panels of Fig.|26]l. Figure|28]shows in detail the 
dihadron azimuthal correlations dNpair/dA(^ in central AuAu collisions II191I |195]| : 
the away-side hemisphere shows a very unconventional angular distribution with a 
"dip" at A(|) « 71 and two neighbouring local maxima at A(|) w 7t ± 1.1 - 1.3. 
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Fig. 28. Normalised azimuthal dihadron distributions, l/N,,-igdNpair/dA(^. Left: STAR data 
in central AuAu (squares) and dAu (circles) for p™'-"^ = 1.3 ~ 1.8 GeV/c and two ranges of 
Pj'^ 11951 . Right: PHENIX results in central AmAk for various p'^'S'"-^™^ ranges [ 191 1. 



Such a "volcano"-like profile has been interpreted as due to the preferential emis- 
sion of energy from the quenched parton at a finite angle with respect to the jet axis. 
This could happen in a purely radiative energy loss scenario due to large-angle ra- 
diation 1 196 1, but more intriguing explanations have been put forward based on the 
dissipation of the lost energy into a collective mode of the medium in the form of a 
wake of lower energy gluons with Mach- ||40l Hi] |42l or Cerenkov-like B2ll43ll44l 
angular emissions. 

In the Mach cone scenario BOl 1411 l42ll . the local maxima in central AuAu 
are caused by the Mach shock of the supersonic recoiling parton traversing the 
medium with a resulting preferential emission of secondary partons from the plasma 
at an angle 9m (Fig. |29] l. Such a mechanism would give access, via Eq. ( fTsT i, to 
the speed sound of the traversed matter. In an expanding plasma, the speed of 
sound changes from q = 1/ (QGP) to Cs ~ ^/O^ (hadron gas) through Cs — 
(mixed phase) in the case of a first-order phase transition. The time-averaged value 
is (cj) = ^ Jq dtCs{t) « 0.3 pTl with a resulting Mach angle 9m = arccos(c5) « 1.3, 
see Eq. ( flSb . in rough agreement with the experimental data. 

In the Cerenkov picture ll42ll43l(44J . it is argued that the combination of the LPM 
gluonstrahlung interference and a medium with a large dielectric constant (n w 2.75 
is needed in Eq. ( fT6b to reproduce the location of the experimental peaks), would 
also result in the emission of QCD Cerenkov radiation with the double-hump struc- 
ture observed in the data. However, at variance with the Mach angle which is constant 
in the fluid, the Cerenkov angle decreases with the momentum of the radiated gluon. 
Such a trend is seemingly in disagreement with the fact that the measured 9^ remains 
relatively constant as a function of pf '"'' lil9T][T99l . In addition, STAR |200| and 
PHENIX 1 201 1 three -particle correlations studies, seem to clearly favour the conical 
over deflected-jets interpretation. 
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Fig. 29. Top: Perturbed energy (left) and momentum (riglit) densities for a gluon moving with 
P = 0.99955 in a perturbative QGP 1 197 1. Bottom: Perturbed energy density (left) and energy 
flux (Pointing vector, right) for a jet with P = 0.75 from an AdS/CFT model 1 198 |. 

Theoretically, the disturbance of the energy-momentum tensor caused by a light- 
quark has been studied in a perturbative plasma II197I as well as for heavy-quark in 
a jy — 4 SYM plasma 02021 11981 . In both cases a clear conical structure as well 
as a strong flow generated along the path of the jet (diffusion wake |40, 203 |) are 
observed (Fig. l29T l. The results are sensitive to the viscosity of the medium. Yet, it 
is unclear if phenomenologically such partonic collective wake(s) and cone survive 
both hadronisation and the final hadronic freeze-out 1 2031 1204] 12051 l206l . Results for 
a pQCD plasma 02O5I indicate that the conical signal does not survive freeze-out: a 
peak at A(|) = 7t appears due to the strong diffusion wake. More involved studies, ac- 
counting for e.g. the plasma expansion and the hadronic phase evolution, are needed 
before a final conclusion can be reached. 

5.3 Pseudo-rapidity correlations: near-side "ridge" 

Figurel30lshows the associated Arj-Ac]) particle yield (down to very low /j^""' > 20 MeV/c) 
for trigger hadrons pj^ > 2.5 GeV/c in pp (PYTHIA simulations) and in central AmAm 
(PHOBOS data) at 200 Ge V. Both distributions show a clear peak at (Arj , A(|)) « (0, 0) 
as expected from jet fragmentation, but the near-side peak in heavy-ion collisions 
features a wide associated yield out to Ar| « 4, referred to as the "ridge" [207 |. The 
existence of such unique long-range rapidity coiTelations in the near-side of the trig- 
ger parton which is, by construction, the least affected by the medium, is puzzling. 
The properties (particle composition, pr slope, intra-particle correlations) of this 
structure are very similar to those of the soft underlying event in the collision 02091 . 
This clearly suggests that the ridge is formed from bulk matter and not from jet frag- 
ments. Though many different interpretations have been put forward (see e.g. 02 101 
for a summary), models that do not require jet triggers for the effect to appear - such 
as e.g."glasma" flux-tubes 021 II or "trivial" modifications of the 2- and 3-particle 
correlations due to radial flow flW] - seem favoured by the data. 
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Fig. 30. Per-trigger associated liadron yield for Pj > 2.5 GeV/c as a function of Ar| and A(|) 
for pp (PYTHIA, left) and 0-30% central AuAu (PHOBOS, right) collisions at 200 GeV (208). 



6 Jet observables in AA collisions 

The measurement of the leading fragments of a jet (single-hadron spectra and/or di- 
hadron azimuthal correlations at high-pj-) in AA collisions has been covered in detail 
in the previous Sections. In this last Chapter, we focus on full jet reconstruction 
in nuclear reactions. The study of the energy and particle-multiplicity distributions 
within a jet issuing from the fragmentation of a quenched parton is a powerful tool 
to study the response of hot and dense QCD matter to fast quark and gluons. 

6.1 Full jet reconstruction in AA collisions 

Experimental reconstruction of jets in hadronic collisions is an involved exercise 112 121 
12131 that requires at least three steps: 

(1) Clustering algorithm: Hadrons belonging to a given jet are measured in the de- 
tector (usually in the cells of hadronic and electromagnetic calorimeters) and are 
clustered together, according to relative "distances" in momentum and/or space, 
following an infrared- and colUnear-safr procedure (see below) that can be also 
appropriately applied to "theoretical" (Monte Carlo) jets. The algorithm needs 
to be fast enough to be run over events with very high multiplicities. Various 
jet-finders exist presently that fulfill all such conditions such as e.g. the kj 112 141 
and SISCone |215j algorithms implemented in the FastJet package II216I . 

(2) Background subtraction: Jets are produced on top of a large "underlying event" 
(UE) of hadrons coming from other (softer) parton-parton collisions in the same 
interaction. At LHC energies, extrapolating from dEj / dT\\^=Q = 0.6 TeV mea- 
sured at RHIC fl47|, one expects a total transverse energy of ^^1 TeV in 1-unit 
rapidity at midrapidity in PbPb. Jet reconstruction is usually carried out with 
small cone radius R = ^/Arp+Acp = 0.3 - 0.5 (or similar A:7'-distances, D) to 
minimise the UE contributions. Indeed, at the LHC in a /? = 0.4 cone one ex- 
pects AEt = nxR^ X 1/ (2%) x dET/dy]\-^=o « 80 GeV, with large fluctuations, 
making it challenging to reconstruct jets below Ej ~ 50 GeV. Various UE sub- 
traction techniques have been developed in combination with the kj M2171|218| , 
UAl-cone 112 191 12201 or iterative-cone 1221 1 algorithms. 

(3) Jet corrections: The energy of the reconstructed and background-subtracted jets 
has to be corrected for various experimental and model-dependent uncertainties 
before comparing it to theoretical predictions. Experimentally, the jet energy- 
scale (JES) is the most important source of systematic uncertainties in the jet 
yield and requires careful data-driven studies (e.g. via dijet and y-,Z-jet Ej- 
balancing in proton-proton collisions). In addition, before a given "parton-level" 
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pQCD calculation can be compared to a measured "hadron-level" jet spectrum, 
one needs to estimate the non-perturbative effects introduced by the underlying- 
event and hadronisation corrections. In pp collisions, this final step is carried out 
usually comparing the results from two Monte Carlo's (e.g. PYTHIA and HER- 
WIG) with different models for the UE multiparton-interactions as well as for the 
hadronisation (string- and cluster- fragmentation respectively). 



(1) Jet clustering algorithms 

In practical terms one usually deals with three types of "jets" (Fig. [31] left)- Exper- 
imentally, a calorimeter jet (aka "CaloJet") is a collection of 4-vectors based on the 
energy deposited in calorimeter towers clustered in pseudorapidity-azimuth accord- 
ing to a given algorithm. [Often nowadays, the experiments use also the momentum 
of the (low pt) charged hadrons measured by the tracking system to reconstruct the 
jet energy with improved resolution, so the name "Calo" is not fully justified.]. At the 
Monte-Carlo generator level, a hadron or particle jet (aka "GenJet") is a collection of 
hadrons issuing from the (non-perturbative) hadronisation of a given parton. Theoret- 
ically, a parton-level jet is what one actually calculates in pQCD. The (non-unique) 
method of linking an initial parton to a set of final-state particles (or other objects 
with four-vector like properties) relies on a procedure known as "jet algorithm". 



SISCone. R=1. f=0.5 | 




Fig. 31. Left: Schema of jet production and measurement [2121 . Right: Reconstructed jets in 
r|-(|) space with the SISCone (top) and kj (bottom) algorithms for a simulated pp event at the 
LHC l222i. 



The goal of a clustering algorithm is to combine hadrons into jets according to 
a given "distance" (radius). Theoretically, such a procedure must be infrared- and 
collinear-safe - i.e. adding a soft gluon or splitting a given parton must not change 
the final number of reconstructed jets. In addition, the jet-finder must not be too 
sensitive to nonperturbative effects - hadronisation, underlying-event (and pileup0 
in pp) - and it must be realistically applicable at detector level (e.g. not too slow). 
There are two broad classes of jet algorithms II2I2I |2T3| : 

Pile-up refers to the additional pp collisions occurring in the same bunch-crossing at high 
luminosity (at J^f = 10^"* cm^^s^' one expects ~25 simultaneous collisions at the LHC). 
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• Cone-type algorithms are "top-down" approaches i.e. they identify energy flow 
into pre-defined cones of a given radius. One sums the momenta of aU parti- 
cles j within a cone of radius R around a seed particle / (often the particle or 
calorimeter tower with the largest transverse momentum) in azimuthal angle (|) 
andpseudorapidityrj, i.e. taking all j such that 



The direction of the resulting sum is then used as a new seed direction, and one it- 
erates the procedure until the direction of the resulting cone is stable. There exist 
various cone jet-finders: JetClu, ILCA/MidPoint, ICone, SISCone, ... which have 
been mainly employed at hadron colliders (see e.g. Fig. [31] top-right). Their main 
advantages are their speed, which makes them easy to implement in triggers, and 
the simplicity of the UE corrections. On the other hand, their particular imple- 
mentations can be messy (seeding, split-merge, "ratcheting", missed or "dark" 
towers, ...) and infrared/collinear safety is not guaranteed in many cases. 

• Sequential clustering algorithms are "bottom-up" approaches that rely on pair- 
wise successive recombinations of the closest hadrons in momentum up to a 
given (predefined) distance D. One introduces distances dij between entities (par- 
ticles, protojets) i and j, and dtB between entity / and the beam (B). The clustering 
proceeds by identifying the smallest of the distances and if it is a dij recombining 
entities / and j, while if it is dtB calling / a jet and removing it from the list. The 
distances are recalculated and the procedure repeated until no entities are left. 
The distance measures for several algorithms are of the form 



the jet-radius parameter (equivalent to R in the cone finders), and p parameterises 
the type of algorithm: kj (p=l) |223 |, Cambridge/Aachen (p = 0) |224|, anti-^j- 
ip = -1) M222II (Fig. [31] bottom-right). On the positive side, these algorithms - 
widely used at LEP and HERA - are explicitly infrared and collinear safe and 
more "realistic" than the cone-based ones as they mimic (backwards) the QCD 
shower branching dynamics. On the other hand, they used to be slow and lead 
to jets with irregular shaped which complicated the UE subtraction compared 
to the cone jet-finders (this has been now solved as discussed in the next subsec- 
tion), making them not competitive in a heavy-ion environment with very large 
hadron multiplicities. Recently, the time taken to cluster particles has been sig- 
nificantly improved in the FastJet |2161 package, based on Voronoi diagrams, 
going down for the default kj jet-finder from G{N^) to 6'{N\nN). Jet clustering 
in nucleus-nucleus collisions is now routinely performed at sub-second times. 

(2) Underlying event subtraction 

Background energy in a jet cone of size R is ff{R^) and background fluctuations are 
&{R). As aforementioned, the soft background from the underlying event in a cone 
ofR = 0.4 in centi-al nucleus-nucleus collisions at RHIC (LHC) is about 40 (80) GeV. 
Fig. [32] (left) shows the (charged) jet and background energies as a function of the 
cone radius R in ALICE |157, 219]. Jets can only be identified if the background 

Yet, with large multiplicities, the kj algorithm (often labelled a "vacuum cleaner") has 
actually an average area --^ nR^, whereas modern versions of the cone finder (assumed to 
have always an area TZF?') with split-merge steps such as SISCone, turn out to be quite non- 
conical, with small areas ~ TtR^ /2 that renders them efficient in noisy environments 12251 . 



A?; = (11,-11;)' + (^,-C^>)'< 



(48) 



where A?- 
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Fig. 32. Left: Charged jet energy in a cone of radius R (full lines) in ALICE compared to the 
background energy from a HIJING I.226J PbPb simulation for different cuts in the particles pj 
(dashed lines) 11571 . Right: STAR AiMm dijet event after background subtraction I227II2281 . 

energy within the cone is smaller than the signal energy. This can be achieved by 
using small cone sizes (e'j^'' oc R^)^ R = 0.3- 0.5, and/or by applying pr or energy 
cuts on the charged hadrons or calorimeter towers. The latter option is not optimal 
since it also introduces potential biases in the investigation of jet-quenching effects. 
STAR Il227l |228| (Fig. [32l right) uses a seeded-cone algorithm with R = 0.4 and 
p™' = 0.1-2 GeV/c, and estimates the UE background event-by-event from the av- 
erage energy in cones without seeds which is then subtracted from the reconstructed 
jets. ALICE uses a modified version of the UAl-cone algorithm (R = 0.4) where the 
mean cell energy from cells outside a jet cone is recalculated after each iteration of 
the cone jet finder and subtracted from all cells 1157 . 219 1. 

Similarly, CMS 11591 122 1 1 subtracts the UE on an event-by-event basis with a 
variant of the iterative "noise/pedestal subtraction" for pp collisions II229II . Initially, 
the mean value and dispersion of the energies in the calorimeter cells are calculated 
for rings of constant pseudorapidity, r|. The value of this pedestal function, P(r|), is 
subtracted from all cells (the cell energy is set to zero in case of negative values) 
and the jets are reconstructed with the default ICone finder. In a second iteration, the 
pedestal function is recalculated using only calorimeter cells outside the area covered 
by jets with Ej > 30 GeV. The cell energies are updated with the new pedestal 
function and the jets are reconstructed again, using the updated calorimeter cells. 

Alternatively, FastJet [217| proposes a background-subtraction procedure af- 
ter running any infrared-safe algorithm. The method is based on the concept of a 
'jet area' A constructed by adding infinitely soft particles ("ghosts") and identifying 
the region in r|,(|) where those ghosts are clustered within each jet |218|. Each recon- 
structed jet pT is then corrected by subtracting the median value of the noise distri- 
bution in the event, p — median [{;?7-/A}], in the jet-area A, via p™'' — pr —A p. 
In practical terms, one fits the measured PT{^) /A background distribution for each 
event with a parabola form, p(r|) — a + br\^ (which excludes any jet peak) and cor- 
rects then the jet pr using the pf'' formula above. 

(3) Jet energy corrections 

The last step of any jet analysis consists in correcting the pr of any measured CaloJet 
to match closely that of the associated GenJet and/or PartonJet, so that it can be 
compared to theoretical expectations. In principle, the different corrections can be 
decomposed as shown in Fig. |33] The experimental coiTections (labelled Level 1-5 
in the plot) can be extracted from the data themselves. For example, the correction 
LI (noise offset) can be obtained from minimum bias events without jet activity, and 
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experimental MC-dependent 




Fig. 33. List of typical factorised jet-energy corrections (CMS analysis) 12291 . 



the L2 (flattening of relative, r] -dependent, pr responses of the calorimeters) and 
L3 (absolute pT calibration) can be derived using pj-balancing techniques in back- 
to-back di-jet and y-, Z-jet events in pp collisions. A precise calibration of the jet 
energy scale (JES) is essential. Given the steep (power law) fall-off of the jet cross 
section as a function of energy and the relatively large binnin^3 of the jet spectra of 
€?(5 - 30 GeV/c), an uncertainty of 10% in the JES can propagate into uncertainties 
as large as 50% (!) in the jet yield at a given pj bin. The L3 correction is thus 
the most important source of experimental uncertainty in any jet measurement. The 
two last connections, L4 (fraction EMF of energy deposited by hadrons in the EM 
calorimeter), and L5 (flavour correction accounting for the different characteristics 
of - and therefore detector responses to - gluon, light-quark and heavy-quark jets), 
can be e.g. obtained in back-to-back y-jet and /7,c-identified dijet events in pp. 

The two "theoretical" corrections (L6-UE and L7-parton) aim at bringing the pj 
of a CaloJet as close as possible to that of its originating parton. They can only be 
obtained from MC simulations that model the effects of final-state-radiation (FSR), 
hadronisation and underlying-event. FSR and hadronisation tend to remove energy 
out of the jet, whereas the UE has the contrary "splash-in" effect. In pp collisions, 
the total shift on a jet pj due to these effects can be approximated by the uncorrected 
sum {dpj) w {5pt)^^^^ + {^Pt)^ J + {^Pt)^^ |230|. The way these effects modify 
the jet energy as a function of the parton pr, flavour and the used cone radius R 
are summarised in Table [3] Whereas the effect of FSR can be in principle computed 
perturbatively, the UE and hadronisation corrections rely on model-dependent de- 
scriptions of multi-parton interactions (MPI) and parton-to-hadron fragmentation. In 
pp collisions, one usually compares the result of PYTHIA and HERWIG - which have 
different MPI and different (string vs. cluster) fragmentation models - to gauge the 
dependence of the measured jet observables on these non-perturbative phenomena. 

Table 3. Main physical effects contributing to a shift (Spr) in the transverse momentum of 
a jet compared to that of its parent parton in pp collisions (cases with '-' do not have any 
dependence at LO) 12301 . 





Dependence of jet (Spj-) shift on 




parton pr 


colour factor 


radius R 


final-state radiation 




c, 




hadronisation 




Ci 


-\/R + ff{R) 


underlying event 






R^/2 + ff{R^) 



In heavy-ion collisions, in-medium FSR and UE are significantly enhanced com- 
pared to pp jets (see cartoon in Fig. |2]i, but at high enough pT the final parton-to- 
hadron non-perturbative fragmentation occurs in the vacuum and should be the same 



The bin-width is not constant in the whole spectrum but given by the absolute pj resolution 
at each bin: relative jet pj resolutions are in the 25%-10% range for PbPb at the LHC U57l 
[T58l[T59l . 
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as in pp. Ideally, the effects of the UE can be controlled embedding MC jets in real 
events, and the influence of hadronisation can be gauged comparing the results of 
e.g. Q-PYTHIA and Q-HERWIG 1 34 1 (adding also eventually medium-induced modi- 
fications of the colour structure of the jet shower evolution). Effects on jet quenching 
observables - which are the ultimate goal of our studies - can then be isolated com- 
paring the results of different parton energy-loss MCs such as e.g. PYQUEN (with 
large out-of-cone elastic energy-loss) and Q-PYTHiA (with its embedded BDMPS 
radiative energy loss). 



6.2 Jet spectra 

The direct comparison of the fully corrected /^r-differential jet spectra in AA and pp 
collisions will provide crucial first-hand information about the nature of the medium 
produced in heavy-ion collisions. The expected ALICE (ATLAS and CMS) jet pj 
range measured in PbPb collisions for nominal integrated luminosities is pj- « 30 - 
200 GeV/c (50 - 500 GeV/c, see Fig.[34]left). A natural generahzation of the nuclear 
modification factor, Eq. (l23T l. for jets ll47l 12311 . 



(n ■ fo--^ - dN^{pT,R^'^\i>f^^)/dydpl 



is a sensitive measure of the nature of the medium-induced energy loss. The steepness 
of the spectra amplifies the observable effects and the varying values of the jet radius 
^max ^jj^ jjjg minimum particle/tower energy sa co™'" will provide, through the 
evolution of /?^(/?7-;/?™™,(M™") at any centrality, experimental access to the QGP 
response to quark and gluon propagation [47 1. If the medium-induced energy loss of 
the parent parton is radiated inside the jet cone, one will observe w 1 at variance 
with the leading-hadron spectra (Raa ^1). On the contrary, important large-angle 
radiation (as expected e.g. in some models of collisional energy loss [231 J ) will result 
in a quenched iet spectrum in heavy-ion collisions (R^^ < !)■ 
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Left: Jet spectra for various PbPb centralities expected at 5.5 TeV in CMS 
= 0.5 nb^') tl59|. Right: Comparison of the vacuum and in-medium jet shapes for 
GeV gluon jets in central PhPb collisions at the LHC 11671 . 



6.3 Jet shapes 

The study of the internal structure of jets - via observables such as jet shapes and jet 
multiplicity distributions -in pp collisions at Tevatron has provided valuable tests of 
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the models for parton branching and soft-gluon emission in the vacuum 02321 . The 
energy degradation of partons traversing a dense QCD plasma will be also directly 
reflected in the modification of such jet observables in heavy-ion collisions ||32ll47l . 
Three variables are useful in this context: 

• the differential jet shape, p(r), defined as the average fraction of the jet pj that 
lies inside an annulus of radius r± 5r/2 (e.g. 5r = 0.1) around the jet axis: 

• the integrated jet shape, l'(r), i.e. the average fraction of the jet pj that lies 
inside a cone of radius r concentric to the jet cone (by definition, ^{r = R) = \): 

• the thrust variable, constructed from the 3-momenta p,- of all jet particles, char- 
acterises the global energy flow: 

^^max„,^Hi:!ld, (53) 

i.e. ^ = 1 (1/2) if the dijet event is pencil-like (spherical), i.e. if all particles 
are aligned (or not at all) along a thrust axis n^^. The projection of all particle 
momenta on the direction n along which the momentum flow is maximal is the 
"thrust major" and on that orthogonal to the plane formed by n^? and n is the 
"thrust minor" 1233 ] . The JEWEL MC predicts a broadening of these (perturbative 
and infrared- safe) =3^aj and distributions inside a dense QGP 1931 . 

Interestingly, from an experimental point of view, all those observables are ro- 
bust against jet energy scale/corrections/resolutions. Medium-modified jet shapes in 
PbPb collisions at LHC energies have been analytically investigated in [32, 47| (see 
e.g. Fig.[34l right). More detailed studies using the recently available jet-quenching 
Monte Carlo's (Section [J!2l l are needed. 

6.4 Medium-modified fragmentation functions 

Due to the coherence and interference of gluon radiation inside a jet (resulting, on av- 
erage, in angular orderin^\ of the sequential branching), not the softest partons but 
those with intermediate energies (£/, ^ E^^^^^'^) multiply most effectively in QCD 
cascades 1234|. This is best seen in the approximately Gaussian distribution in the 
variable t, — \og{Ejg,/ p) — log(l /£) for particles with momentum p in ajet of energy 
Ejf.,, which peaks at the so-called "hump-back plateau" at intermediate ^ « 3 - 4 val- 
ueo (Fig- ED left). In a QCD medium, energy loss shifts parton energy from high-z 
to low-z hadrons and, as a result, leading hadrons are suppressed as seen in Fig. |35] 
(right) where, for increasing q coefficient, the fragmentation function 
is increasingly depleted at high-z. Correspondingly, the number of low-pj- hadrons 
rises, resulting in a higher humped-back in Fig. [35] (left). 



Angular ordering (or coherence) implies 9^,^, ^ ^kip^ ^ ^k^k^ ^ ^ki,k2 ^ where 6*,;,, 
is the emission angle of the primary soft gluon from the direction of the hard parton, Qk^ki 
is that of the softer secondary gluon from the direction of the primary gluon, etc. 
More generally, the peak is at ^ = 0.5 \n{E jg, /Ag^^) at leading order (e.g. ^ = 3 for 
Ejet = 100 GeV) but it is shifted by an amount oc in MLLA 1.48 J . 
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Fig. 35. Left: Single inclusive distribution of hadrons vs. t, = ln{Ejg, / p) for a 17.5-GeV jet 
in e^e^ collisions (TASSO data) compared to MLLA predictions in the vacuum (solid curve) 
and in-medium (dashed curve) 1331 . Right: Medium-modified pion FF for a 100-GeV gluon 
in a medium of length L = 2 - 6 fm (solid-dashed lines) with increasing q parameter L92J . 



Theoretically, the resummed (next-to) Modified Leading Logarithmic Approxi- 
mation (N)MLLA approach describes well, to (next-to)-next-to-leading order 
accuracy, the measured distributions of hadrons %^h{z,Q^) inside a jet (Fig.|35] left) 
down to nonperturbative scales Q^g k, Aqcd ~ 0.2 GeV, provided that each parton is 
mapped locally onto a hadron ("Local Parton-Hadron Duality", LPHD) [235 1 with a 
proportionality factor K w L Various recent promising applications of the (N)MLLA 
approach 133] |35] |36] |37] 12361 have investigated QCD radiation in the presence of a 
medium. 

Photon-jet correlations 

The y-jet (and Z-jet) channel provides a very clean means to determine parton frag- 
mentation functions (FFs) 02371 . In the leading-order QCD processes of photon pro- 
duction (Compton: qg qy, and annihilation: qq gy), because of momentum 
conservation the photon is produced back-to-back to the hard jet, with equal and 
opposite transverse momentum. Since the prompt y is not affected by final-state in- 
teractions, its transverse energy (Ej) can be used as a proxy of the away-side parton 
energy (E^' w E^) before any jet quenching has taken place in the medium. Once the 
quark fragments into a given hadron h, the y-h momentum imbalance variable ll45l . 
z^i^ = —Pj-ft-Pry/IPT-yP, reduces at LO to the fragmentation variable, i.e. z^j^ — Z- The 
FF, defined as the normalised distribution of hadron momenta l/NjetsdN /dz rela- 
tive to that of the parent parton Ej' , can then be constructed using z^j^ or, similarly, 
t, = — ln(z^,J, for all particles with momentum pr associated with a jet in heavy-ion 
collisions. 

ALICE II2381 12401 and CMS II239I have caiTied out simulation studies of the y- 
jet channel, where the isolated y is identified in the EM calorimeter, the away-side 
jet axis (Ai^y^ jet > 3 rad) is reconstructed in the calorimeters, and the momenta of 
hadrons around the jet-axis (/?,>, < 0.5) are measured in the tracker. In the CMS 
acceptance and for E^ > 70 GeV, about 4500 y-jet events are expected according 
to PYTHIA (scaled by the Glauber nuclear overlap) in one PbPb year at the nom- 
inal luminosity. The obtained FFs for photon-jet events - after subtraction of the 
underlying-event tracks - are shown in Fig. |36]for central PbPb. Medium modified 
FFs are measurable with high significance in the ranges z < 0.7 and 0.2 < ^ < 6. 
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Fig. 36. Ratio of medium over vacuum FFs as a function of ^ for quenclied partons obtained 
in Y-jet simulations for central PbPb at 5.5 TeV (0.5 nb"') in ALICE (left) Ii238i and CMS 
(right) 12391 . 



7 Summary 

We have reviewed the main theoretical motivations behind the experimental study of 
parton scattering and jet evolution and fragmentation in the hot and dense QCD mat- 
ter created in high-energy nucleus-nucleus collisions. The phenomenology of parton 
energy loss has been summarised as well as the main experimental results on single 
inclusive spectra and dihadron correlations measured at high transverse momentum, 
mainly in AuAu reactions at RHIC collider energies, up to ^/s^ = 200 GeV. The 
analysis of jet structure modifications in heavy-ion collisions provides quantitative 
"tomographic" information on the thermodynamical and transport properties of the 
strongly interacting medium produced in the reactions. 

Two chapters have discussed in detail two notable experimental results: (i) the ob- 
served factor of ^^5 suppression of high-/??- leading hadrons, and (ii) the strongly dis- 
torted azimuthal distributions of secondary hadrons emitted in the away-side hemi- 
sphere of a high-pr trigger hadron, in central AmAm relative to pp collisions in free 
space. Most of the properties of the observed high-p^ single hadron and dihadron 
suppression (such as its magnitude, light flavour "universality", pr, reaction cen- 
trality, path-length, and y/s^ dependences) are in quantitative agreement with the 
predictions of parton energy loss models in a very dense system of quarks and glu- 
ons. The confrontation of these models to the data permits to derive the initial gluon 
density dN^ /dy w 1400 and transport coefficient q — ^(10 GeV^/fm) of the pro- 
duced medium at RHIC. 

In the last Chapter of this document, we have reviewed the details of jet reconstruc- 
tion in heavy-ion collisions: jet algorithms, underlying-event background subtrac- 
tion, and jet energy corrections. The analysis of jet spectra, jet shapes and the ex- 
traction of medium-modified parton-to-hadron fragmentation functions at low- and 
high- relative hadron momenta promise to fully unravel the mechanisms of parton 
energy loss in QCD matter, in particular at the upcoming LHC energies. The study 
of jet quenching phenomena proves to be an excellent tool to expand our knowledge 
of the dynamics of the strong interaction at extreme conditions of temperature and 
density. 
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